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1.  Introduction 


In  the  neurophysiological  literature,  it  has  been  well  recognized  that 
a  neurcm  cell  is  spatially  extended  and  hence  that  a  realistic  description 
of  neuronal  activity  would  have  to  take  into  account  synaptic  inputs  that 
occur  randomly  in  time  and  at  different  locations  on  the  neuron's  surface 
[(  17  )].  In  this  paper,  we  shall  develop  a  stochastic  model  to  describe 
the  evolution  of  the  membrane  potential  along  the  surface  membrane  of  a  nerve 
cell.  We  will  regard  this  potential  as  a  stochastic  process  or  random  field, 
indexed  both  by  time  t>0  and  by  location  xeX  where  X  will  represent  the  sur¬ 
face  membrane  of  a  neuron.  Our  work  is  based  on  and  extends  the  work  of 
several  authors,  including  Wan  and  Tuckwell,  Riccardi  and  Sacerdote,  J.  Walsh, 
and  G.  Kallianpur  [  21,  18,  20,  14].  In  all  of  these  earlier  works  the 
neuronal  membrane  X  was  represented  by  a  single  point  [  18  and  14  ] 
or  by  an  interval  (0,b]  of  the  real  line.  The  latter  model  -  the  only  spatial 
model  treated  stochastically  so  far  as  we  know  -  has  been  considered  by  Wan  and 
Tuckwell  and  analyzed  probabilistically  by  Walsh.  The  deterministic  back¬ 
ground  for  their  work  is  core-conductor  theory  and  the  one-dimensional  cable 
equation  which  is  adequate  for  situations  which  involve  "longitudinal  dis¬ 
tances  that  are  many  times  the  cylinder  diameter".  Other  choices  of  X  arc 
more  realistic  if  one  is  concerned  with  locations  "close  to  point  sources 
of  current  and  in  problems  where  the  distribution  of  potential  in  a  large 
extracellular  volume  is  of  primary  importance",  I  18,  p.  42,  also  p-Sb). 

In  the  present  work,  the  form  of  X  will  be  quite  general;  for  example, 
it  can  be  any  smooth,  compact,  d-dimensional  manifold.  The  case  d=0  or  1 
includes  the  work  of  the  previous  authors;  d=2  models  the  surface  of  a  neuron 
(useful  also  in  the  study  of  impulses  originating  at  or  near  the  soma)  and 
d«5  is  suitable  for  modelling  the  interior  of  organs  such  as  the  heart. 


1.2 


Before  discussing  the  models  and  results  it  might  be  of  interest  to  give 
a  brief,  though  simplistic,  description  of  the  physiology  of  neuronal  acti¬ 
vity  (see  [17]). 

A  neuron  is  a  cell  whose  principal  function  is  to  transmit  information 
along  its  considerable  length,  which  often  exceeds  one  meter.  "Information" 
is  represented  by  changing  amplitudes  of  electrical  voltage  potentials  across 
the  cell  wall.  A  quiescent  neuron  will  exhibit  a  resting  potential  of  about 
60  mV,  the  inside  more  negative  than  the  outside,  for  reasons  described  below. 
Under  certain  circumstances  the  potential  voltage  in  the  dendritic  tree  will 
rise  above  a  threshold  point  at  which  positive  feedback  causes  a  pulse  of 
up  to  100  mV  to  appear  at  the  base  of  the  dendritic  tree;  this  pulse  is  trans¬ 
mitted  rapidly  along  the  body  and  down  the  axon  of  the  cell  until  it  reaches 
the  so-called  "pre-synaptic  terminals"  at  the  other  end  of  the  neuron.  Hero 
the  pulse  causes  tiny  vesicles  filled  with  chemicals  called  "neurotransmitters" 
to  empty  out  into  the  narrow  gaps  between  the  presynaptic  terminals  and  the 
dendrites  of  other  neurons.  When  these  chemicals  diffuse  across  the  gap 
and  hit  the  neighboring  neurons'  dendrites,  they  may  cause  the  potential 
voltage  in  these  dendrites  to  rise  above  a  threshold  point  and  initiate  another 
pulse. 

The  molecules  which  make  up  neuronal  cell  walls  (and  all  other  living 
membranes)  are  shaped  a  little  like  lollypops  with  two  sticks.  The  candy  end 
is  a  phosphoric  acid  whose  electrical  charge  is  unequally  distributed,  or 
"polar".  This  allows  it  to  mix  with  water  molecules  (which  are  also  polar) 
and  with  ions  such  as  sodium,  potassium,  and  chloride.  The  two  lollypop 
"sticks"  are  long  fatty  acid  chains  which  are  not  polar;  they  dissolve 
readily  in  oil  or  fats,  but  not  in  water.  These  cells  tend  to  arrange  them¬ 
selves  so  that  their  fatty  tails  meet  only  other  fatty  substances  and  only 


their  polar  ends  directly  neet  water  molecules  and  dissolved  ions;  one  quite 
stable  such  formation  is  a  two-molecule-thick  membrane,  with  all  the  fatty 
tails  in  the  middle.  This  two-layer  thick  membrane  of  molecules  which  are 
half  phosphoric  acid  and  half  fatty  lipid  are  called  "phospholipid  bi layers". 

These  bilayers  make  extraordinary  insulators  and,  since  they  are  so  thin, 
extraordinary  capacitors.  If  for  some  reason  a  charge  imbalance  should  arise 
so  that  one  side  of  the  membrane  is  more  positive  and  the  other  side  more 
negative,  the  insulating  property  of  the  layer  would  tend  to  maintain  that 
imbalance  and  the  resulting  voltage  potential. 

In  fact  nearly  all  living  membranes  exhibit  such  a  charge  imbalance  and 
resulting  voltage  potential,  due  to  their  selective  permeability  to  the  pas¬ 
sage  of  various  ions.  Large  protein  molecules  stick  through  the  lipid  bilayer; 
when  several  of  these  come  together  they  can  form  a  passage  through  the  layer 
like  the  staves  of  a  barrel  and  let  one  or  another  species  of  ion  pass  through 
the  membrane.  Different  proteins  have  the  ability  to  pass  different  species 
of  ion  selectively  -  say,  allowing  potassium  to  pass  but  not  sodium  -  and 
this  gives  rise  to  the  electrical  potential.  In  the  cell  at  rest  potassium 
is  about  30  times  more  plentiful  inside  the  cells  than  outside,  while  sodium 
is  about  10  times  more  plentiful  outside  than  in.  This  imbalance  is  maintained 
by  the  ATP  sodium-potassium  pump,  which  acts  like  a  revolving  door.  It  active¬ 
ly  trades  Na  and  K  to  keep  the  potassium  in  and  sodium  out.  Since  the  membrane 
is  more  permeable  to  potassium  than  to  sodium,  the  positively  charged  potassium 
leaks  out  of  the  cell  down  its  concentration  gradient  into  the  surrounding  fluid, 
while  sodium  is  unable  to  leak  into  the  cell;  the  result  is  a  net  negative 
charge  inside  the  cell  of  about  60  mV.  The  voltage  potential  itself  prevents 
more  potassium  from  leaking  out  since  potassium  ions  are  attracted  by  the  rel¬ 
atively  negative  charge  inside  the  cell. 


Returning  to  the  description  of  our  model,  let  ^(t.x)  represent  the  dif- 

i 

ference  between  the  voltage  potential  at  time  t  at  the  location  xeX  and  the 
resting  potential  of  about  -60  mV.  As  time  passes,  5  evolves  due  to  two  sep¬ 
arate  causes: 

(i)  Diffusion  and  leaks:  Depending  on  the  nature  of  X,  the  electrical 

properties  of  the  cell  wall  may  be  approximated  by  postulating  a  contraction 

2 

semigroup  {T^}  on  L  (X,r)  where  F  is  a  suitable  o-finite  measure  on  X.  Tor 
example,  if  X  =  [0,b],  core  conductor  theory  suggests  the  semigroup  correspond¬ 
ing  to  the  diffusion  equation 

II  =  -6C^+  6AC^  (e,6  ^  0) 

with  Neumann  (or  insulating)  boundary  conditions  at  both  ends.  In  neural 
material  like  heart  muscle  in  which  electrical  signals  can  travel  more  easily 
in  some  directions  than  in  others,  the  Laplacian  should  be  replaced  by  a  more 
general  second-order  elliptic  operator. 

(ii)  Random  fluctuations:  Every  now  and  then  a  burst  of  neurotransmittcr 
will  hit  some  place  or  another  on  the  membrane  and  suddenly  the  membrane  poten¬ 
tial  will  jump  up  or  down  by  a  random  amount  at  a  random  time  and  location. 

It  is  believed  that  these  random  jumps  are  quite  small  and  quite  frequent, 
making  it  reasonable  to  hope  that  they  can  be  modelled  by  a  Gaussian  noise 
process;  in  any  case  the  arrivals  at  distant  locations  or  in  disjoint  time 
intervals  are  believed  to  be  approximately  independent,  justifying  their 
modelling  as  a  mixture  of  Poisson  processes  or  as  a  generalized  Poisson  process 
Our  principal  concern  is  to  prove  the  existence  and  uniqueness  of  solu¬ 
tions  to  stochastic  differential  equations  (s.d.e. 's)  that  describe  the  evolu¬ 
tion  of  the  voltage  potential  C  -  a  special  example  of  such  an  equation  is 


(1.1) 


dC^  *  +  5AC^)dt  ♦  dX^ 

(B,6iO)  in  which  the  "noise"  or  driving  term,  is  a  generalized  Poisson  process 
on  X  or  possibly  a  Gaussian  process  on  R^x  X  -  and  to  prove  that  under 
comparatively  mild  conditions  a  sequence  of  solutions  to  this  equation  with 
Poisson  driving  terms  will  converge  in  distribution  (in  the  sense  of  weak  con¬ 
vergence  of  the  induced  probability  measures)  to  the  solution  with  a  Gaussian 
driving  term. 

These  results  are  established  in  Sections  2  and  3  of  this  paper.  In  our 
formulation  of  the  problem,  the  voltage  potential  C  is  viewed  as  a  stochastic 
process  taking  values  in  4',  the  dual  of  a  suitably  chosen  countably  Hilbcrtian 
nuclear  space.  (Thus  the  stochastic  differential  equations  we  consider  govern 
nuclear  space  valued  processes.)  It  turns  out  that  almost  all  the  paths  of 
the  voltage  potential  process  lie  in  the  Skorokhod  space  D(3R^:  H  ^)  where 
H  ^  is  a  suitable  Hilbert  space. 

The  work  of  Wan  and  Tuckwell  as  well  as  its  rigorous  treatment  in  a  more 
general  set-up  by  Walsh  is  discussed  in  some  detail  as  Example  2  of  Section  4. 
It  is  perhaps  appropriate  here  to  remark  briefly  on  the  relationship  of 
Walsh's  work  to  ours.  As  in  Wan  and  Tuckwell,  Walsh  takes  X  to  be  an  inter¬ 
val  [0,bl  and  considers  the  potential  as  a  stochastic  process  of  the  two  para¬ 
meters  t,  the  time  and  x,  the  location.  The  techniques  of  2-parameter  martin¬ 
gale  theory  are  used.  The  approach  adopted  in  our  paper  leading  to  stochastic 

differential  equations  in  infinite  dimensional  spaces  (specifically  Hilbert 
spaces  or  nuclear  spaces)  has  the  advantage  that  the  theory  can  also  be  applied 
to  more  general  cases  in  which  multiparameter  martingale  methods  are  either 
cumbersome  (e.g.  if  X  £F*^,  d>2)  or  inapplicable  (e.g.  if  X  is  a  sphere  or 


a  more  general  compact  smooth  manifold  with  or  without  a  boundary) .  Ex¬ 
ample  3  of  Section  4  and  Section  5  are  devoted  to  these  applications. 

Though  our  methods  (in  contrast  to  those  of  Walsh)  do  not  permit  us  to 
study  the  behavior  of  the  membrane  potential  process  at  individual  points 
(t,x)  we  are  able  to  prove  stronger  approximation  results  (see  Section  4). 

Stochastic  differential  equations  in  infinite  dimensional  spaces  have 
been  intensively  studied  in  recent  years  in  the  context  of  many  physical 
applications.  They  occur  in  the  work  of  Dawson,  Miyahara,  Holley  and  Stroock 
[  9,  16  and  10].  In  the  last  named  work,  which  is  a  study  of  infinite  parti¬ 
cle  systems  in  statistical  mechanics  the  authors  are  led  to  an  s.d.e.  driven 

d  ^ 

by  )-valued  Brownian  notion.  Similar  equations  were  discussed  by  Ito 

in  his  Evanston  lecture  some  years  ago  [12]. 

The  models  discussed  in  the  present  paper  do  not  investigate  the  pheno¬ 
menon  of  excitability,  or  high-speed  passage  of  voltage  pulses  called  "action 
potentials".  (Xir  models,  however,  are  useful  in  describing  the  sub-threshold 
behavior  of  a  membrane  potential.  Of  particular  interest  is  the  probability 
distribution  of  the  length  of  time  until  the  base  of  the  dendritic  tree  reachc 
a  critical  value  6,  at  which  time,  an  action  potential  is  generated.  Informa¬ 
tion  on  this  is  given  by  our  weak  convergence  result: (Theorem  3.2)  which  im- 
Pl  ies  that  the  first  passage  times  for  the  Poisson-driven  process  converge 
in  distribution  to  the  first  passage  times  for  a  generalized  Ornstein-Uhlen- 
beck  process. 

In  Section  6  we  mention  more  realistic  models  of  neuronal  activity  which 
introduce  nonlinear  semigroups  and  s.d.e. 's  of  ♦'-valued  processes  with  cor¬ 
responding  weak  convergence  theorems.  These  questions  as  well  as  the  problem 
of  estimation  of  parameters  of  interest  will  be  considered  in  a  later  paper. 


2.  Mathematical  Formulation 

Let  ®  self-adjoint  contraction  semigroup  on  the  Hilbert  space 

2 

H  =  L  (X,8,r)  for  some  o-finite  positive  measure  space  (X,8,n.  X  is  intended 
to  represent  a  mass  of  excitable  tissue  (perhaps  the  cell  membrane  of  a  neuron), 
(T^)  the  evolution  semigroup  describing  the  decay  of  the  difference  between 
the  actual  voltage  potential  at  the  time  t^O  and  the  resting  potential  on 
X.  The  measure  T  has  no  physical  significance;  it  is  chosen  for  convenience  in 
order  to  make  {T^}  self-adjoint  and  to  satisfy  several  assumptions  below. 

We  shall  require  that  the  resolvent  =  Jp  e  compact  (for  each 

a>0’) ,  and  even  that  it  satisfy  the  assumption 

^1 

A1 :  For  some  rj>0  the  operator  (R^)  is  Hilbert-Schmidt. 

By  the  Hille-Yosida  theorem  {T^}  has  a  negative-definite  infinitesimal  generator 
-L.  Since  R^  =  (a+L)"^,  the  assumption  implies  that  H  is  separable  and  that  L 
admits  a  complete  orthononnal  set  of  eigenvectors  in  H  with  eigenvalues  0<Xj^5X2S.. 
satisfying 

-2r 

r(a+X^)  ^  . 

These  properties  hold  for  any  a>0  if  and  only  if  they  hold  for  all  a>0r.  we  shall 
take  (1=1  and  set 

(2.1)  0j  *  i:(i+x^) 


Denote  by  <•,*>  the  inner  product  on  H  and  let 

=  {(^cH;  E  <  (1+Xj)^^  <  “»  for  all  r  f  F  } .  For  each  real  number  r 

define  a  quadratic  form  <  and  norm  11*11^,  on  $  by 


<  »  Z  <  <|(,((»^>  >  (l+Xj)^’^ 

Mj.  =  (< 


(2.2) 


a’ld  let  be  the  Hilbert-space  completion  of  $  in  the  inner  product  <  *,*  >^. 

Give  4>  the  Fr^chet  topology  determined  by  the  family  {11*11^}  of  norms  and  let 
be  with  the  inductive  limit  topology. 

The  following  are  straightforward  consequences  of  A1  and  the  definitions 
above: 

i)  $  is  contained  in  the  domain  of  l"  for  every  integer  n, 

ii)  L  and  map  $  into  itself. 

iii)  Finite  linear  combinations  of  {4*^}  are  dense  in  $  and  in  every 
are  orthogonal  in  each  H^. 

iv)  Hq  =  H. 

v)  ^  H  3  H  3  $  if  -«  <  r  <  s  <  “;  the  injection  of  H  into  H  is 
Hilbert -Schmidt  if  s  2  r  +  r^^. 

vij  H  ^  and  are  in  diiality  under  the  pairing  for 

^  c  H  (t>  e  H^.  The  pairing  is  symmetric  and  independent  of  r,  and  (|>T  = 
for  C,(j)  e  H. 

vii)  may  be  identified  with  the  dual  space  (in  the  weak  topology)  to  4>. 

-tX. 


The  proof  of  i) ,  for  example,  proceeds  by  showing  that  the  sum  Z<4),(t'j>e 


4*; 


-1 


converges  in  $  for  each  t20  and  each  e  H  to  T  d)  and  that  t  ((1)-T  <j>)  = 

.1  -‘N 

>  t  (1-e  is  Cauchy  in  H  =  Hg  as  t-+0,  with  limit  L(j>  =  Xj(^^ 

Statements  ii)-  vijaxe  easy,  and  vi.)  is  a  consequence  of  the  nuclear  theorem  (sec 
GV,  Ch  1  Thm  3). 

Let  p  be  a  a-finite  measure  on  RxX  satisfying  the  assumption 


A2:  The  bilinear  form 

0(^,i|))  *  /  a^<ji(x)i|;(x)p(dadx) 
on  $  X  $  is  continuoxis 

and  let  mc^'.  By  the  nuclear  theorem  there  exist  numbers  r2  e  R  and 
0,(  R  such  that 

(2.3)  +  Q(<)),<|))  s  e-ll(j)|l^  for  all  4)  c  4»  . 

i  r2 

All  processes  and  random  variables  are  assumed  to  be  defined  on  a  fixed,  but 
arbitrary  complete  probability  space  (fl,F,P).  Let  N  be  a 


2.3 


I 


regular  Poisson  measure  on  IR  X  x  with  mean/covariance  measure  p(dadx)dt 
and  define  a  ^'-valued  stochastic  process  with  stationary  independent  incre¬ 
ments  hy 

(2.4)  X  r(J)]  =  tm^c]))  +  /  a())(x)l'N(dadxds)-p(dadx)ds'>  . 

^  FxXx(0,tl 

The  interpretation  is  that  N(AxBx(0,t])  should  be  the  number  of  voltage  pulses 

of  sizes  aeAcR  arriving  at  sites  xeB^X  at  times  s<t;  the  probability  that 

-X  k 

exactly  k  such  pulses  arrive  during  the  indicated  period  is  e  X  /k!  with 
X  =  ty(AxB).  A  computation  will  verify  that 

iX  r<l>]  •  1 

(2.5)  logEe  =itmr(j)]+t/  -  1  -  ia<|>(x))p(dadx) 

for  all  t  6  ,  (j)  e  4>. 

Now  let  V  =  V  +  f  e  be  the  voltage  potential  on  X  at  time  t^O, 

L  K  L 

V_e4>'  the  resting  potential,  and  C*  =  V  -  V_  their  difference:  we  will  model 
as  the  ^'-valued  solution  with  initial  value  -  Vj^  to  the  stochastic 

differential  equation 

(2.6)  d?^  =  -L'C^.dt  +  dX^  . 

Here  and  below  L'  and  denote  the  adjoints  of  L  and  when  regarded  as 
operators  on  In  Theorem  2.1  below  we  construct  a  solution  to  (2.6)  by 
evaluating  the  stochastic  integral 

(2.7) 

Before  doing  so,  we  pause  to  introduce  three  examples;  they  will  be  developed  i 

more  detail  in  section  4  below. 

fixamples: 

1.  X  has  a  single  point. 

In  this  case  H  =  Rand  T^  is  multiplication  by  e~^^  for  a  constant  L>U 
Obviously  A]  is  satisfied  for  any  r..  >  -®. 


2.  X  is  an  interval,  L = -A  +  6  (with  Neumann  boundary  conditions). 


Here  we  may  take  Lebesgue  measure  for  T  and  verify  A1  for  all  >  1/4. 

3 

3.  X  is  the  unit  sphere  in  R  ,  L = -A  (Laplace-Beltrami) . 

If  we  take  T  to  be  surface  measure  then  A1  is  satisfied  for  all  Tj  > 
since  then 

-2r-  “  -2r 

0  =Z(1+X.)  =  I  (2£-H)  (!+£(£+!))  ^  . 

^  ^  1=0 

In  all  three  examples,  condition  A2  is  satisfied  for  any  measure  p  ‘’f  the 


(2.8) 


y(AxB)  =  ?  1  (abv^CB)  +  I  1  (-af)vf(B) 

k=l  ^  ®  ®  £=1  ^  ^  1 


in  which  {a^}  c  (0,“)  are  the  possible  sizes  of  ’’excitatory"  (i.e.  positive) 

pulses,  {-a^}  the  sizes  of  "inhibitory"  (i.e.  negative)  pulses,  and  {v^.Vjl  arc 

finite  measures  on  X  giving  the  local  arrival  rates.  The  proof  hinges  on  the 

fact  that  in  each  case  the  orthonormal  set  {<[).}  satisfies  a  bound  of  the  form 
r  ^ 

|(t)^(x)|  s  c(l+Xj)  uniformly  in  x  e  X  and  j  <  “;  we  may  take  r  =  0  in  examples  1 

and  2,  and  r  = -s  in  example  3. 

k  £  k  ^ 

Rxample  1  (with  and  point  masses  of  size  f^  and  fr,  respectively  on 

the  one-point  space  X)  appeared  in  [TO  and  fGK)  and,  with  slightly  different 

notation,  in  Trs].  Example  2  (with  X  =  rO,b]  and  y  satisfying 
2 

■^RxX  ^  y(dadx)  <  «>  appears  in  FJW],  as  a  generalization  of  the  example  in  rwT  !. 

k  £ 

In  the  latter  example  y  was  of  form  2.8  with  each  and  a  point  mass  at 
the  point  x^.  Example  3  was  suggested  to  us  as  a  model  for  excitation  at  the 
soma  of  a  neuron.*  The  basis  for  this  example  is  the  fact  that  synaptic  inp\its 
may  occur  also  in  the  somatic  region  coupled  with  the  usual  assumption  about  th 
approximate  spherical  shape  of  the  soma.  We  know  of  no  examples  other  than 
these  three  which  have  been  studied  to  date.  We  will  motivate  and  introduce 

*bv  Dr.  T.  McKenna  (personal  communication) 


classes  of  new  examples  in  Section  5. 


Consider  now  the  problem  of  performing  the  stochastic  integration  indicated 
in  (2.7)  for  a  specified  random  element  of 

Since  any  element  (|)  c  il>  may  be  expanded  into  a  series  of  the  form  ^ 
which  converges  in  each  (and  therefore  in  "I*),  we  may  hope  to  write  - 

I  <  (J),())^  >  with  given  by 


(2.9) 


Kl-e 


t  -(t-s)X. 
f  e  J  dx^r^,.  ] 


In  order  to  carry  out  this  program,  introduce  the  notation 

^0  ""  ’ 


(2.10) 


m-'  =  mC4).)  , 
J 


xj  =  . 

yj  *  xj  -  tm^  , 


Mj  =  / 


t  sX. 

/  e  ^  dvj  . 
0  ^ 


1  1  k 

Note  that  Y-^  is  a  square-integrable  martingale  with  covariance  function  lEY-'Y  = 

sX.  ^  ^ 

min(s,t)Q((t)j,<l)p  (because  has  independent  increments)  and  e  ^  is  trivially 
predictable,  so  the  martingale  integral  in  (2.10)  is  well-defined  and  is  also 
a  square-integrable  martingale.  We  can  and  do  take  versions  of  and  y]  (and 
therefore  of  X^)  with  sample  paths  in  the  Skorokhod  space  D(R^:1R)  of  right-con¬ 
tinuous  F- valued  functions  on  F^  with  left-hand  limits  at  every  point  of  (0,°"). 


Lenunn  2.1 


-tX.  . 


For  each  T  >0  and  q^r.  +r2,  E  sup  (e  ^M^)^(l+X.)~^^  ^  16T0,f 

^  ^  3«»  0<t5T  ^  ^  ^ 

In  particular,  the  sum  converges  almost  surely. 


Proof: 


so  we  have  the  bound 


(t-s)X 


s  2  sup  Ivjl. 

OssSt  ^ 

Dooh's  inequality  applied  to  yields 

Esup  (e  ^mJ)  s  4  E  sup  |yJ| 
t<T  ^  t<T  ^ 


s  16E  (y|)^ 


=  16  T  J  a^(|)j  (x)vi(dadx) 
2r- 

s  16  T  e2(i+x^)  ^  . 

It  follows  by  the  monotone  convergence  theorem  that 

E^  sup  (e  + 

j  tsT  ■’ 


ZCr.-q) 

S  16T02i:(l  +  X^)  ^ 

s  i6Teje2  , 


and  hence  the  series  converges  almost  surely. 


Theorem  2.1 


For  each  random  element  Cq  of  satisfying 
A3t  " 


for  some  r^  €  R  ,  the  series 


(2.12) 


c,  = 


converges  uniformly  in  0<tsT  in  the  H  topology  for  each  T>0  and  q^inax(r,  +  r-,-T 

-q  1  i 

to  a  process  C,  whose  sample  paths  lie  in  the  Skorokhod  space  D(F^:  H  ^)  of  right - 
continuous  H  ^-valued  functions  on  with  left  limits  at  each  point  of  (0,“'). 

The  process  satisfies 


and  (a  fortiori) 


for  some  C.^  <  ». 


Esup  lU  r  ^  C 
t<T  ^  “H  * 


JE  |C^r<|)]|2  s  Cj.(|(})||2  (OitST) 


Proof:  For  each  T>0  and  n<n'  c  tJ  the  triangle  inequality  yields 

(Initial  term) 


sup  II  I  ^  sup  II  j;  e  j  cj-l).! 

0<t<T  n<jsn'  J  T  n<t<T  n+i  J  M 


Ost<T  n+1 


O^tsT  n+1 
n 

I 

OstsT  n+1 

Wc  bound  the  three  terms  separately: 


n'  -tX.  . 

+  sup  11  ^  e  ^  mJ4)  IL 

n<r<T  T14.1  *-  J  M 


sup  II  I  mj(j)  IL 

1<+<T  n*l  *•  J  M 


(Martinqale  term) 


(Mean  term) 


-tX. 


(Initial  term)^  s  sup  I  (e  ^  C  ^(1  +  X .)  ■^‘’ 
OstsT  j>n  0  J 


4  ^  J 


0  a.s.  (as  n-**)  since  E  r  <  • 

0  -q 


(Martingale  term)^  5  sup  ^  (e  ^  M^)^(l  ♦ 

OstsT  j>n  ^  ^ 

-►  0  a.s.  as  n-*-®  by  Lenuna  2.1. 

(Mean  term)^  s  sup  I  C/q  e  ^  m^ds) ^(1  ♦  X.)'^^* 

Ost^T  j>n  ^ 

j>n  •*  ^ 

^  *7 

-►  0  as  n-^®  since  ||b  H  s  ||ib||  <  0,  <  ®. 

-q  -T2  ^ 

Thus  the  partial  sums  converge  uniformly  on  ro,T]  in  the  H  ^  topology 

Since  each  partial  sian  lies  in  D  =  D(F^  :  H  this  implies  that  the  partial 
sums  form  a  Cauchy  sequence  in  the  complete  metric  space  D;  let  C,  denote  their 
limit . 

The  estimates  above  and  the  Cauchy-Schwartz  inequality  together  yield 

E  sup  lU^ll  i  3E  r(Initial  term)  ♦(Martingale  term)^  +  (Mean  term)  1 

OstsT  ^ 

5  3[03  +  16T0J02  ♦T^02] 


Theorem  2.2 

For  each  OSr^t<®  the  process  constructed  in  Theorem  2.1  satisfies  the 
equation 


=  T;  C  ♦  /,  *,T*  dX 

t  t-r  r  ^ (r,t)  t-s  s 

satisfies  for  each  ()>  e  $  the  equation 


1  ,e. 


(2.131 


-2.9- 


Proof: 

For  any  q2max(r.  ♦r-,  -r_)  the  following  series  converge  for  all  <l>  ^  by 

X  M 

Theorem  2.1: 


tX^  sX.  .  -tX.  .  -(t-s)X. 

^  dY^  e  ^  ^  ds)  , 


c,r*i  .  r  ♦r^.jKcJ  e 

-(t-rlX.  .  -rX.  sX.  -rX.  -(r-sV 

=  Z  e  J{Cje  JdY^e  ^  e 

.  -tX.  sX.  .  -tX.  .  -(t-s)X. 

=  E  (lirtl*  J  (cj  e  ^+/(0,r]®  ^■""^■^0®  ^  ds)  , 


hence 


-{t-s)X.  .  -(t-s)X.  . 


^  m^ds) 


■(t-s)X 


)  jvJ 
s 


f,  dX  PT^  ())1  . 
“'(r.t]  s  t-s’^ 


n 


Corollary  2.1 

Lot  =  oCCq.X^:  Ossst}  be  the  smallest  a-algebra  over  which  and 

X  r^tl  are  measurable  for  all  sst  and  <^c  $.  Then  for  rst, 
s 

E  rc^r<^]!F^]  =  ’"’’'^t-s*''  ' 


Corollary  2.2 

has  the  (strict)  Markov  property,  i.e.  f  ^  is  conditionally  independent 
of0{^^r(|»1:  s2r,  <|>c$}  given  a{ 5^(4))  :  4>e4>}. 

Proof : 

[,et  ’Pjj  e  4>  and  2  r  for  isksK.  Using  the  independent  increments  property 
of  X,.it  is  possible  to  compute  explicitly 

iE^t  t 

logEPe  ^  I  ^  ^  r'^k^  */r'^'"'^’^t  -s^  ^ 

^  kSK  ^  K  r  s 

*  f  re^^nx,s) 


-  1  -  iaF(x,s)p(dadx)ds  , 


where  FCx.s)  =  I  J  since  this  is  not  only  IF -measurable  but 

kSK  ’  k'  * 

o{Cj.r4'l:  e  -measurable,  we  are  done.  11 

In  the  remainder  of  this  section  we  state  and  prove  results  similar  to 
Theorems  2.1  and  2.2  for  stochastic  differential  equations  of  the  form 

dC^  * 

in  which  the  $' -valued  stochastic  driving  process  is  Gaussian.  The  proofs 
are  similar  to  (and  a  bit  simpler  than)  those  for  the  generalized  Poisson  pro¬ 
cess  already  considered,  so  they  will  only  be  skeletal. 

Several  authors  have  considered  infinite-dimensional  Wiener  processes  and 
Omstein-Uhlenbeck  processes  such  as  W^  and  in  case  ^  =  S.  For  example,  Tto 
has  an  excellent  account  in  his  1981  Evanston  lecture,  '‘12].  Miyahara  considers 
similar  processes  in  connection  with  a  vibrating  string  problem  ri61;  Holley 
and  Stroock,  in  their  treatment  of  a  problem  involving  infinite  particle  systems 
also  introduce  these  processes  flO].  However,  since  no  discussion  in  the  litera 
ture  seems  to  include  all  of  the  estimates  we  will  need  in  Section  we  prefer 
to  derive  them  here. 

Now  let  Q  denote  a  continuous  positive-definite  bilinear  form  on  (j>,  let 
m€<^',  and  let  be  a  path-continuous  ^'-valued  independent -increment  stochas¬ 
tic  process  with  characteristic  functional 

Tto  r  12)  calls  such  a  process  a  "Wiener  S'-process"  when  $  =  S.  ITie  Minlos  and 
nuclear  theorems  allow  us  to  construct  such  a  process  with  continuous  sample- 
paths  lying  in  H  ^  for  any  qSr^  +r2  if  m  and  Q  satisfy 

A2':  |mr<|il|^  +  Q(<^,<|>)  s  e_|I(^||^  for  all  (|>  r  <I)  . 


For  random  ^  H  set 


nd  i;  1  =  e  +  e  •'  m;  +  ml  . 

t  0  t  t 

Theorem  2.3 

For  each  random  element  of  4>'  satisfying 

A3*:  ®  r  “  ®3  r^  e  H  , 

3 

the  series 

(2.15)  = 

converges  uniformly  in  OstsT  in  the  H  ^  topology  for  each  T>0  and  q2max(rj+r 
to  a  process  whose  sample  paths  lie  in  the  space  C(F^  :  ^  q)  of  continuous 
n  -valued  functions  on  F  .  The  process  satisfies 

-q  + 

E  sup  IK.  11^  S  C 
OstST 

and  (a  fortiori) 

E  K^r(D)|^  s  C.J.1I(|>|1^  (Ost^T) 
c_  =  sre.  +  i6Te.8.  ♦T^e.i  <  »  . 


for 


Proof: 


In  this  Gaussian  case  each  z]  is  an  ordinary  Wiener  process,  with  diffusion 
rate  E(Zj)^  «  it  follows  that  Z^,  w|,  and  aay  all  be  taken  to  have 

continuous  saaple  paths.  This  forces  the  partial  subs  I  to  have  continu- 
ous  saaple  paths.  The  proof  that  the  sequence  of  partial  suns  converj^es  uni- 
fomly  on  each  fO.T]  proceeds  just  as  in  Theoren  2.1,  using  exactly  the  same 
estiaates  and  leading  to  the  saae  bound  Cj.*  F! 

Theorea  2.4 

For  each  0Srst«»  the  process  constructed  in  Theorea  2.3  satisfies  the 
equation  t]"'^t-s  ****^s'  for  each  ^  e  <t>  the  equa¬ 

tion 

(2.16)  ♦  /y 


Corollary  2.3 

Let  » 

)  -  y'^Vr^’  *  * 

Corollary  2.4 

has  the  strict  Markov  property. 

Proof : 

The  theorea  and  its  corollaries  are  proved  in  exactly  the  saae  manner  as 
Theorem  2.2  and  its  corollaries.  (1 

Remark  2.1 

If  r,„  is  jointly  Gaussian  with  {W^}  then  C,  will  be  a  Gaussian  H  ^-valued 
process  which  is  in  many  ways  an  infinite-dimensional  Omstein-Uhlenbeck  (hence¬ 
forth:  0-U)  process.  If  m»0  and  Cq  is  independent  of  W  then  for  each  j 


is  the  usual  0-U  process  satisfying  the  stochastic  differential  equation 


(2.17) 


with  initial  value  Cq* 


For  X.  >0  the  process  will  be  stationary  if  is  independent  of 

3  V  u  * 

and  has  a  normal  distribution  with  mean  0,  variance  Q(()).  ,(j).)/2X.;  if  X.  =  ()  then 

3  3  3  3 

®  Wiener  process  starting  at  . 


Remark  2.2 

The  processes  will  be  independent  if  the  are  independent  and  in¬ 
dependent  of  {W^r(j)];  taO,  <!>€$}  and  if  diagonalizes  Q,  i.e.  if  Q(il).  ,(f>j)  =0 
for  i^j;  independence  may  fail  in  other  cases.  Similarly,  the  processes 
constructed  in  Theorem  2.1  will  be  orthogonal  if  ,(|>.)  = /a^(l)j^(x)(|). (x)vi(dadx)  =0 
for  i/j  (as  will  happen  when  r(dx)  =  /a^Vi(dadx),  for  example)  and  if 
E  IY^ri|)3  »  0  for  all  taO  and  €  $. 

Remark  2.3 


The  mean  functional  m  in  the  definition  of  (similarly  J;^)  plays  an  inessen¬ 
tial  role;  in  most  cases  it  can  be  absorbed  into  Cq  and  Vj^  as  follows. 

Assumption  A1  implies  that  L  has  a  finite-dimensional  null-space  spanned  by 
{(j)^ , . . . ,  for  some  n  2  0;  let 

"•o  =  I  in[<|».J(}). 
jsn  ^  ^ 


be  the  projection  of  m  onto  that  null  space,  and  set  m 

K 

converges  in  H,  and  satisfies 
1-r, 


(mVx.)^..  The  series 
j>n  ^  ^ 


m  *  mQ  ♦  ' 

We  can  now  write  m^  (see  (2.11)  and  (2.14))  in  the  form 


if  j  s  n 


.  -tX. 

[m^Cl-e  ^)/Xj  if  j  >  n 

-tX. 

injjr(|i  Jt  +  Cl  -  e 


and  rewrite  ®  as  Vj^  +  tm^  +  C^  where  +  is  the  new  resting  poten 

tial  and  -  nij^  -  tuip  is  the  new  deviation.  The  process  with  initial 

«*«#  #s» 

value  ^0  “  Vq  -  Vj^  =  satisfies  the  equation 


(2.18)  ft':'*’  • 

This  is  just  (2.7)  with  ^  replaced  by  and  by  the  process  -  tm 

which  satisfies  an  equation  like  (2.4),  but  with  m  =  0. 

The  effect  of  m  is  to  bring  the  resting  potential  Vj^  to  a  new  equilibrium 
V|^.  Usually  L  has  no  zero  eigenvalue,  so  does  not  appear. 

Remark  2.4 

Theorems  2.2  and  2.4  suggest  that  the  equations 
<2^  =  -L’C^dt  +  dX^ 
and 

d^t  =  -L'C^dt  + 

might  have  stationary  solutions  on  -oo<t«*»  satisfying  (2.13)  and  (2.16)  for 
all  -«'<r<t«».  When  X^  >  0  for  all  j  this  is  indeed  the  case.  After  setting 
m  =  0  (see  Remark  2.3)  and  defining  ®  o*'  \  ®  tc  F  ,  it  is  easy 

to  construct 

^*(<l>)  =  dx  rr 

;_oo  s  t-S^ 

or 

■  /!» 

as  in  Theorems  2.1  and  2.3. 


.  Weak  Convergence  of  Solutions 


Now  let  us  fix  a  o-finite  measure  space  (X,D,r),  a  self-adjoint  contraction 
2 

semigroup  {T^}  on  H  =  L  (X,8,n,  and  an  initial  voltage  distribution  and 

consider  the  effect  of  varying  the  distribution  of  the  incident  noise  process 
X^.  This  distribution  is  uniquely  determined  by  the  mean  functional  m'^  <t>’  and 
the  measure  y  which  gives  the  expected  frequency  with  which  impulses  of  various 
magnitudes  hit  X  at  various  points. 

In  Theorem  2.2  we  have  derived  the  stochastic  differential  equation  satis¬ 
fied  by  the  "electrotonic  potential"  (=  the  difference  between  the  voltage 
or  membrane  potential  and  its  resting  value  .  As  was  shown  in  Theorem 
2.1,  C,  is  a  process  whose  sample  paths  lie  in  D(R^  :  q) •  To  study  the  con¬ 
tinuous  approximation  of  a  sequence  of  such  processes  it  is  first  necessary 
to  derive  auxiliary  results  on  the  weak  convergence  of  D(R^  :  H  -valued  pro¬ 
cesses.  We  begin  with 


Lemma  3.1 

Let  {P**}  be  any  family  of  Borel  probability  measures  on  D  =  D(F^:H)  for 
some  real  separable  Hilbert  space  H  with  norm  11*11.  Then  {P*^}  is  tight  if  and 
only  if  for  each  e  >  0  and  T  >  0  there  exist  a  finite-dimensional  subspace  F  c  H 
and  positive  numbers  b,6  such  that  each  P°  satisfies 

(3.1)  i)  p“{h€D:  sup  II  (I -njh^ll  >  e)  <e  , 

0<t<T  ^  ^ 

ii)  P®{h  c  D:  sup  Ijllp  h  ||  >  b}  <  e  , 

OSt^T  ** 

iii)  P^^ChcD:  sup  ||n  (h  -h»)l|>  e)  <  e  , 

0Ssi6  ®  " 

P°‘{heD:  sup  1(11  (h- -  h  )  H  >  e)  <  e  , 

T-fisssT  ® 


and 


iv)  P™{h  f  D:  sup  min(  ||llp  (h^  -  h^  )  II  ,  PpCh^  -  h^)  II  )  >  e}  <  r  . 


The  supremum  in  iv)  extends  over  all  choices  of  (tj,t,t2)  satisfying 
and  t2-tjS6.  In  i)  -  iv) ,  IIp  is  the  orthogonal  projection  in  H  with  range  F. 

Proof : 

By  Prohorov's  theorem  it  suffices  to  prove  that  a  set  KcD  has  compact  clo¬ 
sure  if  for  each  e  >  0  and  T  >  0  there  exist  b  >  0,  6  >  0,  and  F  <- H  as  above  such 
that  every  h  e  K  satisfies 

i)  sup  II  (I  -n  )h  11  <  e  , 

Ost<T 

ii)  sup  ||n  h  II  <  b  , 

0<tsT 

iii)  sup  ||n  (h  -  h„)  II  <  e,  sup  l|np  (h^- h  )  1|  <  e  . 

iv)  sup  min{||np  (h^  -  h^^)  H  ,  pp  (h^  -  h^)  11 }  ^  e  • 

The  space  D  with  its  Skorokhod  topology  is  a  complete  metric  space  under  the 
metric 

d(h,k)  =  sup  e"^minCl,  inf  max( |||  Xl|L,  sup  ||h  -  k.  ||)1 
0<T«»  XeA  0St<T  ^ 

where  A  is  the  space  of  strictly  increasing  maps  from  to  and,  for  T>(), 

and  XeA,  l|lX|ll„  =  sup  llog[(Xt  -  Xs)/(t  -  s)3l .  The  proof  is  similar  to  that 
*  05s<t^T 

of  Theorem  14,2  of  Billingsley  rb]  .  The  useful  property  of  this  particular 
metric  is  that  for  every  T  >  0  it  satisfies 

(3,2a)  d(h,k)  <  max(e“^,  dClj-^  ,p^h,  IpQ 

and 

(3,2b)  d(h,k)  ^  max(e“^,  sup  l|h  -  k  |I)  . 

0st<T  ^  ^ 

For  each  n  e  IN  let  and  6^  be  the  subspace  and  positive  numbers  pro¬ 

mised  in  (3,1)  for  e  =  e~"  and  T  =  n;  without  loss  of  generality  we  may  take 

r  c  F  b  sb  ,,  and  6^6  ,  for  each  n.  This  insures  that 

n  n+1  ’  n  n+1  n  n+1 

-  >1 


(3.3) 


llUp  hll  ^  l|n  hll  for  n  <  m 


for  each  h  e  H. 


Let  {h^}  c  K  be  any  sequence  and  define  n^^  =  i.  By  Theorem  15.2  of  ''6  1 

there  is  a  subsequence  along  which  lip  h^  converges  in  D(rO,ll:  F)  and  hence 

along  which  lip  Ij-^  j-|h^  is  Cauchy  in  (D,d).  For  each  j  ^2  there  is  a  further  sub 

sequence  n..  along  which  lip  Ij.^  .^h^  is  Cauchy;  set  =  h^  and  consider  the  di- 

]i  ^  l  ,3  n  ^  1  n^^ 

agonal  subsequence  {k^}.  We  claim  {k^^}  is  Cauchy  in  D. 

-N 

Fix  any  n  >  0  and  take  N  e  K  large  enough  that  e  s  n/3.  Let  F  =  and  fix 
N*  >  N  so  that  for  all  integers  i,jiN*, 

From  (3.2)  it  follows  that 

(3.4)  d(n_k';'.  n„k';')  <  max(e'^.  n/3)  =  n/3  . 

r  1  r  3 

Since  every  n^^  with  i  2  N  is  of  the  form  nj^^  for  some  j  si,  (3.1)-(3.4)  yield 
for  each  i ,  j  s  N* 

d(k!,  k^!)<d(k!,  n_k^)  +  d(II_k^  n_k^.)  +  d(n^k^ ,  k^!)  s  n  . 

'  i’  3^  '  1  F  1  F  i’  F  j  F  j’  3'^ 

Since  each  sequence  {h^} c  K  has  a  Cauchy  subsequence,  K  has  compact  closure.  H 

Lemma  3.2 

Let  {P™}  be  a  family  of  probability  measures  on  D  =  D(11^:  H)  for  some  real 
separable  Hilbert  space  H  with  norm  ll'H.  A  sufficient  condition  for  {P°^l  to 
be  tight  is  that  for  every  T  >  0  and  c  >  0  there  exist  a  finite-dimensional  sub- 
spacc  FcH  and  positive  numbers  c^,  c^,  Cj  such  that  for  every  0<tj^t^t2-T, 
each  P^  satisfies 

i)  J  sup  II  (I-n  )h  ll^dp"  5  e, 

^  0<ssT  ® 

ii)  /  sup  i|n  h  ll^dP®  <  c.  , 

^  0<ssT  ^  ®  ^ 


'in  /nl|np(h  -h  )||^  S  c,(t--t  )  , 


iv)  |p||np(h^-h^^)||^  ||np(h^^-hp||^dp“  <  C3(t2-tj)^  . 

lii 

Proof : 

Chebyshev's  inequality  yields  i)-iii)  of  Lemma  3.1  while  the  same  argument 
Billingsley  uses  in  the  proof  of  his  Theorem  15.6  r6]  gives  iv) .  II 


Lemma  3.3 

Let  Gj  and  G2  be  subsets  of  a  topological  vector-space  with  compact  closures 
Gj  and  G2.  Then  ^  has  compact  closure  as  well. 

Proof : 

The  set  G^  ■*•^2“  *^2‘  ^i^^i’  is  the  continuous  image  (under  ad¬ 

dition)  of  the  compact  set  G^  G2<=VxV,  and  so  is  compact.  This  closed  set 
contains  G^  and  hence  G^,  which  must  therefore  be  compact.  M 

Corollary  3.1 

If  the  families  {PoX"^}  and  {P  «Y"^}  of  Boxel  measures  on  V  induced  by 

a  01 

random  elements  {X^:aeA}  and  {Y^  :  a  e  A}  are  both  tight,  then  the  family 

{p  o  z"^  }  induced  by  Z  =X  +Y  is  also  tight  (even  if  X^  and  Y  are  not  inde- 
a  '  a  a  cx  a  a 

pendent) . 

Proof : 

Fix  e  >  0  and  find  compact  sets  G^  c  v  satisfying  IP  (X^  e  G^ ] > 1  -  e/2, 

P  TY  e  G  1 2 1  -  e/2  for  all  a  e  A;  then  P  rz  €  G  ) > 1  -  e,  where  (as  above) 
a  2  a  j 

G3=  (gj  +g2'-  " 

Theorem  3.1 

Let  A  be  any  index  set  and  let  {m“  :a€A}c4)»  be  an  equicontinuous  family 
of  linear  functionals  on  (u'^-.acA}  a  family  of  measures  on  P  x  X  such  that 
the  bilinear  forms  q“((|),i|j)  =  J  a  <t)(x)i|;(x)y‘^(dadx)  are  equicontinuous  and,  for 
each  ({)€$,/  a^(J>^(x)u“(dadx)  is  bounded  independently  of  a  e  A  by  some  number 
e  ((b)  <*.  Let  be  a  random  element  of  satisfying 


for  some  >  -«>.  Then  for  all  sufficiently  large  qe  IR  the  family  :  ae  A} 
of  Borel  probability  measures  induced  on  D  = D(F^  ;  H  by  the  processes 
{r,  :  a e  a)  constructed  in  Theorem  2.1  is  tight. 


Proof: 

By  the  equi continuity  condition  there  exist  numbers  62  >  0  and  r2  e  R  such 

that 

(3.5)  lm“((l))l^+Q“((|).(t.)  seJl(|)ll^ 

2 

for  all  4)  e  4'  and  ot  €  A.  Let  q  >  max(rj  +  r2,  -r^)  and  set  6^  =  E  llCQll_q  <  “’•  Tt 
follows  that  assumptions  A2  and  A3  are  satisfied  for  each  aeA,  so  Theorem  2.1 
establishes  the  convergence  of  the  sum 


in  the  D  topology  (and,  in  fact,  uniformly  in  the  H  norm  on  O^tiT  for  each 
T  >  0) .  By  Corollary  3.1  it  suffices  to  show  that  the  families  {P?;  •’i=l,2,i; 
induced  by  the  sums 


-tX. 


i  =  l:  ^  ^0 


-tX. 


i  =  2;  =  Le  ^  . 

t  ^3  t  ^  j 


i  =  3:  m 


a 


arc  each  tight. 

i  =1:  Since  does  not  vary  with  a,  the  "family"  {Pj}  consists  of  a  single 
inner-regular  Borel  measure  on  D,  and  hence  is  tight. 


i  =  2:  Fix  e  >  0  and  T  >  0.  By  Lemma  3.2  it  is  enough  to  find  F  and  c. .  c  ,  c,  inde 

J  ^  O 

pendent  of  a  satisfying 

E  sup  11(1  -  n  <  e 


i) 


3.6 


ii)  E  sup  ^  c 

0<s<T  s  1 

iii)  ‘'21^2-hl 

iv)  B||np(MjX,>ll?qll’'F'«?2-<)ll!q  *  's' *2  ' ‘l ' '  ' 

By  Lemma  2.1  we  can  satisfy  i)  and  ii)  by  setting  c.  =160.02  and  letting  F  ho 

-2r^ 

the  space  spanned  by  {(j), , . . . ,  <|).,}  for  any  J  large  enough  that  16T0^y.  ^(1  +A.)  < 

1  J  J 

Now  fix  any  05tjSt2^T  and  compute 

!|n  (M^  )ll^-  =11  I  (e  wj H ?« 

^2  ’^l  ^  j<J  ^2  ^1  J 


-(t--s)X.  .  ^  -(t,-s)X.  •  „  , 

=  I  (I+XJ'^^^C  J  e  ^  dY^»“-  J  e  ^dY^’  ) 

jSJ  ^  (0.t2]  (0,t,] 


(O.tjl 


=  I  j  f  (s)dY^»“)^ 

jSJ  J  (0,«)  J 

-(t2-s)X.  -(tj-s)X. 

for  the  function  f.(s)  =e  -e  1^^^^^., 

Note  that 


(s)  . 


(3.6) 


lf.(s)l.X.|t2-tJl(0^^^3(s).l^^^^^^^(s)  . 


Straightforward  computations  with  the  log  characteristic  functional 

iY^r  4>) 

log  Ee  ^  -  1  -  ia4»(x))v”(dadx) 

2  4 

show  that  for  f,g  in  L  (E^)  nL  (E^)  and  4>,\|>  in 

E  /  f(t)dY“r<f.]  =  0 


(3.7) 


E(/  f{t)dY“r(|,])2  ^  I  a^(),^(x)f^(t)p"(dadx)dt 


Q°‘(4>,<I>)J  f^(t)dt 


(3.8)  E(/  f(t)dY^r4,])2f/  g(t)dY“ri|;])^  =  (/  a‘*/(x)i|;^(x)p“(dadx)/  r(t)e/(t)dt) 


.Otr  .  2 


4^2,  ,.,2, 


2  2 


5. 


+  2(/  a^<|)(x)iKx)pVadx))^(/f(t)g(t)dt)^ 

+  (/  a^(|>^(x)u“(dadx)J  f^(t)dt)(|  a\^(x)vi“(dadx) J  g“(t1df) 
=  /  a^(|)^(x)4)^(x)p“(dadx)/  f^(t)g2(t)dt 
+  2(Q°‘(<|,.,p))2(/  f(t)g(t)dt)^ 

+  Q“(4>.<l>)Q“(<li. '!>)(/  f^dt)(;  g^dt) 
s  /  a'‘<l)^(x)t|)^Cx)M“cdadx)/Cfg)^dt 
♦  3q“c<I>,(J))Q“(<|).iP)(/  f^dt)C/  g^dt)  . 

For  the  particular  case  of  4i  =  (J)^  and  ^  =  f j  w®  find  (by  (3.6)  and  (3.7)) 
E(/  f.(t)dYj*“)^  =  f?(t)dt 

^  Q°‘(<j...(j).)rtjX2|t2-tj|2  +  |t2-tj|] 

2r 

se2(i+x^)  +  n|t2  -  tjl 

so  iii)  is  satisfied  with 


-2q+2r2  2  2 

c  =0  I(i  +  x.)  ^rrxSn 

jsj  J  J 


-(t-s)X  .(t  -s)X 

If  we  put  (|)  =  4)^,  ip  =  4>^  ,  fj(s)  =  e  (s)  -  e  ^(0  t  ) 

-(t2-s)X  -(t-s)X 

«j(s)=e  ^^(0,t]^*^"®  ^  IjQ  (s),  then  by  (3.6)  and  (3.8), 


lEllnp  IlHp ir 


=  I  (1 +X  )-2‘^(14X f  (s)dY^*“)2(J  g  (s)dYj*'^2 
i,j5j  J  s 

^  I  (1  +  X  )-^‘’(l  +  X ^"^‘’(JaV^Jdy^^dfJg^dt) 
i,jsj  J  ^  J  ^  J 

43  I  (1  ♦X 4X  )'%“(0  .4.  )0“((j)..4).)(/f^dt)(/g"dt) 
i.i<J  ■'  ^  ^  J  J  1  .1 


3.8 


<  In* \.y\i  *  ^  xjt  )  (t  - 1.)  ct,  - 1) 

ijsj  ^  J  1  J  ^  J  Z 

X  T 

+  3  I  Cl*X.y\l*\  y^%hl*\  )  ^(l*\)  2(x2t2^1)(t-t  )fX^t5M)(t,-t) 
i,j<J  J  *  ^  1.1.- 


"  C3lt2-tjl 


=  V(  I  e^((t).))rAJ+TxJj(i +  V4[eje2(i +1^^)  V 


Thus  i)-iv)  are  satisfied  and  {P™}  is  tight. 


ot  ot 

i =  3;  Since  is  nonrandom,  each  is  concentrated  on  a  single  point.  It  is 
sufficient  (by  Lemma  3.1)  to  find,  for  each  e>0  and  T>(),  a  positive  number  (S 


such  that 


i)  sup  II  (I  -n  )m“||  s  e 
0<tST  r  c  q 

ii)  sup  !|n  m^ll  <00 

iii)  sup  ||n  rm“-m“]l|  <  e 

0SSS6  F  s  0  -q 

sup  ||n  rnL-m“]|l  5  e 

T-dssST  F  ‘  s  q 

•  _  _  _ •  _  y  n  TT  r___Q^  tII  11' 


iv)  sup  min(llnprm^-m^^]ll_^,||nprm^^-m^l|l  )  <e  . 

By  (3.5)  we  have  {m'^}cH  uniformly  bounded  by  ||m'^||^  ^0,;  pick  .1  r  W 

-r2  -r2  ^ 

2(r2~q)  2  2 

large  enough  that  (1  + X^)  5  e  /(02T  )  for  j  > J  and  let  F c  H  ^  be  the  space 

generated  by  {(|)j^,...,  <))j}.  It  follows  that 

sup  11(1-  n  )m“l|^  =  sup  I  (1 +  X.)"^‘^(m“r(t  1/  e  ds)^ 

OSt^T  r  e  -q  Q^tST  j>J  •’  0 


S  T^(gW0,)  I  (1 +X.)  ^m“r(|)  1)- 

j>J  J  ^ 


.  (eVe2)ll«  II!, 


verifying  i) 


i  e 

The  same  estimates  show  that 

sup  lln  "“11^  5  T^llm""!! 

OstsT  ^  ^  ^ 


s  T 


2 II  a  I 


'  -r 

*  *1 


2 


<  <» 


so  iL*)  holds. 

For  iii)  and  iv)  it  is  sufficient  to  prove  that  Hpm^  is  a  uniformly  continu 
ous  function  of  t  on  each  finite  interval,  i.e.  to  find  6 > 0  so  that 

(3.9)  llnprm“-m®]ll_^s  e  if  OSsstsT,  |s-t|s6. 

Wc  compute 

lln  rm“-m“]||2  =  I  (l^X  )’2‘?(m“r<|>.]J  e  du)^ 

^  ®  jsj  J  J  s 

^|t-s|2||m«l|2^ 
s82|t-sl^  , 

so  (.3.9)  holds  for  all  T>0  with  6®e/(02)**.  □ 

It  follows  from  Theorem  3.1  that  {P°*}  has  at  least  one  cluster  point,  and 
that  any  cluster  point  is  a  Borel  probability  measure  on  D.  We  now  turn  to  a 
situation  in  which  there  is  a  unique  cluster  point  to  which  every  sequence  in 
{P*}  mu.st  converge. 


Proposition  3.1 

Let  Y  and  Z  be  independent -increment  -valued  processes  with  log  charac¬ 


teristic  functionals 


3.10 


(3.10)  log  Ee  ^  =  t/  -  1  -  ia(|)(x))u(dadx) 

iz  r(j)] 

log  Ee  ^  =  -}stQ(<|),<|)) 

for  some  O-finite  measure  y  on  R  X  such  that  the  positive-definite  bilinear 
forms  Q  and 

=  !  a^<j>(x)if»(x)y(dadx) 

satisfy  the  bounds 

(3.11)  Q(4..()))  ^  0||(j)||J 

Q^(<l>.<l>)  i  elkllj 

for  some  numbers  0  e  and  re  R  .  Extend  Q  and  to  by  continuity. 

2 

Then  there  are  unique  continuous  linear  maps  Y  and  Z  from  L  (R^  ;  H^)  to 
the  square-integrable  random  variables  such  that 

(3.12)  Y^r4,]  =  YC1^q^^3<|.)  . 

• 

These  maps  have  log  characteristic  functionals 

iYrpl  f 

(3.13)  '{'^(F)  =logEe^’'‘  *^-'  =  J(e  ®  -  1  -  iaF^ (x))y(dadx)ds 

^^(F)  =logEe^^'’’’^  =  -hj  Q(F^,F^)ds 
which  satisfy  the  inequality 

(3.14)  lY(F)  -  'J'(G)  I  <  Js0(||1f|1|^  +  lllGlIlp  |||F  -  G|1|^ 

^  O  1 

for  F,Ge  L^(R^  :  H  p.  Hero  1|1f||1^  denotes  the  norm  (J  l|Fj|^ds)^  of  an  element 
Ft  L^(R^  :  H^). 

Proof: 

For  functions  FeL^(R  :  H  )  of  the  form  F  (x)  =  ^  f.(s)(J>.(x)  with  each  f. 

+  r  s  3  3  3 

a  step  function  with  compact  support,  (3.13)  follows  immediately  from  (3.10). 

By  taking  derivatives  in  (3.13)  one  sees 


BYrp]  *  Ezrp]  >  0  , 

F.(Yrp32)  *  Q^(P^,P^)ds  <  e|||F|||J  . 

F(zrp]^)  =  QCFs.Vds  S  e|||Pll|2  . 

2 

Since  step  functions  P  of  the  fom  indicated  are  dense  in  L  (F^  :  H  p,  the 
proof  will  be  complete  (by  continuity)  once  we  verify  (3,14)  for  step  functions. 
First  we  consider  Y.  If  x  and  y  are  any  two  real  numbers  then 

I -  1  -  ix)  -  (e^^  -  1  -  iy) I  =  ( /yjj  i^e^®  dsdt 1 

^  IJ^Itldtl 
^  ’sdxj  ♦  |y|)|x-y|; 

hence 

IH'y(F)  -yG)l  S  %/a^(lFl  ♦  |G|)lF-Gldvds 

s  J*(/a^(|F|  +  |G|)^dlJds)’*(/a^|F-G|^dvds)’* 

^%0|ll  |fMg|  IIIJIIf-gIII^ 

^  ’se(|||F|||^.|l|G|||^)|||F-GlIl^  . 

Pot  Z  we  use  the  parallelogram  law: 

-^^(G)!  S  Js/lQ(F^.Fp  .Q(G^.G^)lds 

»  Js/lQCF^+G^,  F^  -G^)|ds 

sMjQCFg  +  G^,  F^+G^)ds)**(jQ(F^-G^,  T^-G^)ds)^ 
^  >se||lF+G|||^|||F.G|I|^ 

^  %0(II|f|||^HI|g|||^)|||f-g|I|^  .  □ 

Remark  3.1 

For  each  tj»e  and  t  >  0  the  random  variable  appearing  in  the  proof  of 

■Ihcorem  2.1  (resp.,  2.3)  is  of  the  form  Y^F)  (resp.,  zrp])  with 


^  (j).  (x)  .  It  is  easy  to  verify  that  Fe  I.  (F  :  H  ) 
j  +  r 


if  since 

<  ty.4)r<().]^(i  + 

^3  3  3 

=  tlkll^ . 


Remark  3.2 

For  each  \p  e  t  and  t  >  0  the  number  m^rip]  appearing  in  the  proofs  of  Theorems 
2.1  and  2,3  satisfies  the  inequality 


-(t-s)X. 

-  t’’ll"II.JI|F|ll^  . 

For  each  nc  K  let  x"  be  an  independent  increment  ^'-valued  process  with 
log  characteristic  functional 

ixjrcj)]  . 

log  Ee  =  itm"r<|)]  +  tj(e^®*^'-  ^  -  1  -  ia(J>(x))y’^(dadx) 

for  some  m*' c  $'  and  o-finite  measure  y”  on  Rx  X  for  which  the  bilinear  functional 

=  /  a^(t)(x)ip(x)y"(dadx) 

is  continuous  on  Let  w”  be  an  independent-increment  ^'-valued  process 

with  log  characteristic  functional 


log  E  e  *  itmr(}»]  -  istQ((j),(l>) 

for  some  mc$'  and  continuous  positive  bilinear  form  Q. 

Let  {c!!)  and  all  be  random  elements  of  'l>'  satisfying  the  bounds 


(J.ltl) 


I 


3.13 


biu;ii^,^®3 
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for  some  fixed  Tj ^  and  6^ €  (independent  of  n) ,  and  such  that  each  Cy  is 


independent  of  {x”}  and  Cy  is  independent  of  {W^}.  Let  and  ^  be  the  4>'- 
vulued  processes  given  by 


(3.17) 


=  T'E”  ♦  f  T'  dx’’ 
^  t^O  ■'(0,tVt-s  s 


*  /(o.ti’-i-s  “"s 


in  Theorems  2.1  and  2.3,  respectively.  We  now  prove  the  main  result  of  this 
section. 


Theorem  3.2 


In  addition  to  (3.16)  assume  that  (p”).  On"},  and  (Cq)  satisfy  the  following 


conditions: 


A4; 


i)  Ee 

ii)  mr<t)l 


i^_r<j)]  idrd)] 

=  lim  Ee  " 


n-H» 


=  lim  m"r<|>} 

n-Ko 


iii)  q(<|),(p)  =  lim  Q"(<t>,i|)) 

n-Hx> 


iv)  0 


=  lim  /  a^(j)^(x)v"(dadx) 


n-»oo 


V)  (m"r4).}|^  +  Q"(<(),<|))  S  e 


for  every  (j)  and  ip  in  4>,  where  independent  of  n.  Then  for  every  q  > 


max(rj  *^2,  -r^)  the  measures  p"  =  Po(5")~^  induced  on  D  =  D(F^:H  by 


-1 


converge  weakly  to  P  *  P°(C) 


Proof: 


From  ii)  and  iii)  it  follows  that  m  and  Q  satisfy  v)  too,  i.e.  that 

|mr(J)]|^  +  Q(4).<l>)  ^  0j|4>ll^ 

for  all  (J)€$.  By  Lebesgue's  dominated  convergence  theorem  it  follows  that 

(3.18)  |lm-m"ll2^  =  I- (•»(<)>.]- m"r,j,.])2(i  + 

-*•0  as  n-*-®  , 

(3.19)  llQ-  Q"ll?q«.q  =  (Q(<t)i  .4>j)  -  Q"(«J)i  )  ^(1  +  “^^(1  +  X  .) 

-►0  as  n  . 

The  norm  ||s||^g  =  S((|).  ,<|>j)^||4»j  ||J||(j)^  ||^)'^  appearing  in  (3.19)  is  the  tensor 

product  norm  1|*||^  ®  ||*llg  of  the  bilinear  form  S,  when  S  is  regarded  as  an  element 
of  (4i«4>) '  =  ^r*^s  [IQj.lef.  43.3  and  the  discussion  following). 

and  satisfies  |S((|),(|))I  S  IIs|| . 

Thus  the  hypotheses  of  Theorem  3.1  are  satisfied,  so  {p”)  is  tight.  To  prove 
that  p"  converges  weakly  to  P  it  suffices  to  show  that  the  finite-dimensional 
distributions  of  p”  converge  to  those  of  P,  for  then  {p")  could  have  no  cluster 
point  other  than  P  and  so  must  converge. 

Fix  any  K  e  IN  ,  {tj, . . . ,  tj^}  c  ,  and  {ipj . c:  <I>.  It  remains  only  to 

show  that 

(3.20)  Ee  =  lim  Ee 

n-w 

2 

Fix  any  T  s  max{tj^}  and  define  F  e  L  (F^:  H^)  by 

note  that 

Ilf,!!,  ^  Jkll^’kl'r 


for  any  reF  and  Fg  =  0  fo*"  s  i  ro,T],  so  Fq  ^  and  (||f||(^  s 

As  in  the  proofs  of  Theorems  2.1  and  2.3  we  decompose  X,  and  C*'  into  the 


three  series 


(t-s)X. 


-tX.  -(t-s)X. 

each  series  converging  uniformly  (for  tsT)  in  the  H  ^  topology.  Thus 


-(t-s)X 

X  J  , 


■{t-s)X 


-t,  X. 
It  j 


(t,  -s)X. 
k  1 


-(t,-s)X. 


=  *  ‘^ro,T]®^^s^‘*® 


and  similarly 


Jk^V^kl  =  *  Iro.T]  • 

k 

2 

where  Z  and  y”  are  the  continuous  linear  maps  from  L  (F^:  H^)  into  the  square 
integrable  random  variables  given  in  Proposition  3.1. 


By  independence  we  have 


(3.21) 


'i'  =  log]Ee  "  =  '{'^^(Fq)  +  'F^(F)  +  fJF) 


'i'r  (PJ  = 
^0 


'{'^(F)  =  logEe 


izrp] 


=  -Wrn  TiQCPe.FJds  . 


Similarly, 


(3.22) 


f  =  logEe 


=  n""’  ‘  ’’  n'”’ 

^  1  m 


'*'  =  log  E  e 

^0 


iaF  (X) 

'!'„(?)  =  J(e  ^  -  1  -  iaF  (x))y  (dadx)ds 

s 

nf^)  ■  ‘/ra.i/r^W^  • 

m 

Fix  G  >  0  and  choose  J  e  IN  large  enough  that  the  orthogonal  projection  IT  of  H 
onto  the  span  of  <|)j}  satisfies 

IJICI-  n)i(-Jl^  s  e/(1202T^lllF|||^)  ; 

by  the  triangle  inequality  and  Remark  3.1  we  have 


(3.23)  |!|(I-n)F|||^  ^  ’^^'f'kllq 

^  G/(1202||lF||lq). 

By  i) ,  iv) ,  (3.18)  and  (3.19),  we  can  find  Ne  N  large  enough  that  for  each  n>N 


{3.24a) 


(3.24b) 


(3.24c) 


(3.24d) 


llogEe  ^  ^  -  logEe  ^  ^  |  <  e/3  , 

||m-  m"||_q  <  e/3T*^lllFll|^  , 


II Q-  Q^'lLqm-q  < 


e/12ll|F|ll‘  , 


I  /aY(x)y"(dadx)  <  e^/4T0  |||f||1^(  J,  (  I  l!)'^’)^ 

j^J  J  ^  q  j<j  k<K  ^ 


3.17 


Introduce  the  temporary  notation  G  =  IIF,  i.e. 

-(t,-s)X. 

G  (x)  =  I  I  ^  1^  -,Cs)4).(x).  Then  ||1g|||  s  |||  f|||  . 

®  j<J  k<K  ^  ^  ^ 

and 

SO  /aV(xly”(dadx)ds  <  T(L(Lli;/.  I  /a'^d)  Cx)%"(dadx) 

^  s  3  K  K  3  jsJ  J 

<  e^/402lllFl(l^  by  (3.24d). 

Thus 

jf-CG)  -  T  (G)  I  =  |/(e^^^-l-iaG+?5a^G^)dv"ds-*sJa^G^dy"ds+*s/Q(Gg,G^)ds| 

£* 

5  /le^^^-l-iaG+JsaV|dy"ds+Js/|Q"(Gg,Gp-Q(Gg,Gg)  |ds 
5  ijlaGlWds  +  JllQ''-  QlLq*.q  IIgJI^  ds 
s  |(/aVdy"ds)^(JaVdy"ds)*^+  1|q"-  Qll.q^.^^  ll|G||lq 

<  e/12  +  e/12 
=  e/6  . 

By  Proposition  3.1  and  (3.23), 

(3.25)  If^CF)  -  nCn  I  ^  I’l'^CF)  -  'V^iG)  \  *  |'i'2(G)  -  T^j^(G)  1  +  |'i'^^(G)  -  ^^^(F)  1 

^  ’'S02(lllF|llqHl|Gll|q)l||F-G|||q+e/6+%02(|||F|||q.|||G|||^)|||F-Gy 

<  e/6  +  262l||F|||^|||F-nF|||^ 


<  e/6  +  e/6 
=  e/3. 


3.18 


By  (3.21),  (3.22)  and  (3.24b)  we  have 

(3.26)  -  \(n|  5  JlmrF^] -Tn^rPgllds 

m 

<  /II 

<  llm-  m"!)  T^'^dllEjl^ds)’'' 

=  llm-  m"|l_^T''1||F|l|^ 

<  e/3  . 

Finally  by  (3.21),  (3.22),  (3.24a),  (3.25)  and  (3.26), 

|¥-'1'"|  <  |f^  (Fq)  -'J'  ^(Fq)|  +  If^CF)  -H'  j^(F)|  +  |'{'^(F)  -T  ^(F)  | 

^0  Cq  Y  m 

<  e/3  +  e/3  +  e/3 
=  e  . 

This  completes  the  verification  of  (3.20)  and  the  proof  of  Theorem  3.2.  II 

Remark  3.3 

Condition  A4  iv)  can  be  replaced  by 

A4  iv)  '  lim  /  la  |  ^*'^  |(|)(x)  |  ^^'^p"(dadx)  =  0 

n-xw 

for  some  5>0.  Since  in  most  applications  the  incident  impulses  arc  uniformly 
bounded  by  a  constant  A < «>  and  A4(iv)'  and  A4(iv)  are  equivalent  in  this  case, 
we  omit  the  necessary  changes  in  the  proof  of  Theorem  3.1. 


4.  Examples  and  Applications 

In  applications,  the  set  X  represents  a  neuron  or  some  part  of  a  neuron 

(such  as  the  cell  wall,  soma  or  dendritic  tree)  or  an  assemblage  of  neurons 

and  muscles  such  as  the  heart.  In  this  section  we  discuss  several  possible 

mathematical  models  for  X,  the  associated  measure  F  and  the  self-adjoint  con- 

2 

traction  semigroup  on  L  (X,r)  which  models  the  decay  of  voltage  potentials 
on  X  in  the  absence  of  arriving  excitatory  and  inhibitory  impulses.  We  will 
also  consider  classes  of  impulse  arrival-rate  measures  {y^}  for  which  we 
can  verify  the  hypotheses  of  Theorems  2.1,  3.1  and  3.2  and  so  construct  solu- 
tions  C  to  the  equation 

dC"  =  -L'C"  dt  +  dx" 

(see  Section  2)  and  verify  that  converges  weakly  to  a  Gaussian  process  f,. 

First  we  consider  in  greater  detail  the  three  examples  introduced  in 
Section  2.  It  is  convenient  to  state  and  prove  the  following  useful  lemma. 

2 

Lemma  4. 1  Let  be  a  semigroup  on  H  =  L  (X,r)  satisfying  A1  and  also 
(4.J)  c,  =  sup  sup  |<|>  (x)|  (1+A.)”’"  <  <» 

X  j  ^  ^ 

for  some  r  <  «>.  Then  any  o-finite  measure  y  onlRxX  satisfying 

(4.2)  c  =  /  a^y(dadx)  <  «> 

FxX 

also  satisfies  assumption  A2. 

Proof;  From  (4.1)  and  (4.2)  it  follows  that 


satisfies 


=  fa^(p(x)>p(x)ii(dadx) 

^  /a^Cjd+A.)’^  Cj(l+Xj^)’'y(dadx) 


and  hence 


Example  1.  X  consists  of  a  single  point. 

Without  loss  of  generality  we  may  take  r(X)  =  1  and  identify  each  ij)  c  <J> 

with  its  value  <|i(X)  €  IR,  in  other  words,  we  identify  ^  with  F.  Any  contraction 

semigroup  is  of  the  form  T^({)  =  e  (J)  for  some  L  e  IR^  ,  and  each  is  one- 

dimensional  with  ||<|)||g  =  |<}>(X)|  (1+L)^.  Each  measure  y  on  IR^X  is  determined 

uniquely  by  its  marginal  y.(A)  =  /  y(dadx)  on  the  Borel  sets  in  IR,  and  sat- 

^  A 

2 

isfies  A2  if  and  only  if  /^a  y^Cda)  <  “•  For  measures  of  the  form(2.8)(with 
p  <  00  and  q  ^  ,  we  have 

(4.3)  Pj(A)  I,(-a^') 

where  f^  =  v^(X)  and  ff  =  v?'(X) ,  aind  A2  becomes 
e  e  1  1 


(4.4) 


This  is  trivially  satisfied  when  p  and  q  are  both  finite. 


A  sequence  p  of  measures  on  IR*X  satisfies  conditions  (iii)-(iv)  of  A4 
Lf  and  only  if,  for  some  number  i  0,  the  marginals  y^(*)  =  p^(*xX)  satisfy 

(4.5)  /a^p'J(da) 

/a%j(da)  -*■  0 

as  n  for  of  form  (2.8)  a  sufficient  requirement  is 


,,  r  Jc,n,  k,n.2  r  r^,n,  i,n.2  2 

(4.6)  ^  f;*  (a^*  )  >  (a. '  )  -  o 

max{a*^’”,  -►  0 

k.il  ®  ^ 

as  n  In  either  case,  for  any  sequence  {m"}  c  IR  there  is  a  sequence  of 

unique  (in  distribution)  generalized  Poisson  processes  x”  =  tm*^  +  fax 

n  n  ^  IRx(0,t]  ^ 

(N  (dads)-Pj(da)ds)  as  in  (2.4),  satisfying  (2.5)  and,  for  each  sequence  Cq 

of  square -integrable  random  variables,  a  unique  (by  Theorem  2.2)  solution  c” 

to  the  stochastic  initial  value  problem 


d?"  =  -LCj  dt  +  dX^ 


c"  =  c" 


If  p"  is  of  the  fonii  (2.8)  then  X*)  may  be  represented  as 


(4.7)  X"  =  t(m"-lii")  >  I  a;'"N;’"(t)-  I  a^V’"(t) 


k,n  „k,n, 


i,n„«.,n, 


k=l 


)l=l 


where  N^’"(k,S.=  l,2, . . .)  are  independent  Poisson  processes  with  variance 

6  1 


parameters  f^’*^  and  f^'*'  and 


(4.8)  in"  =  ^  fX,n^k,n  _  ^ 

I,  ®  ®  a  i  i 


assuming  the  latter  sum  converges;  for  this  it  is  sufficient  that  each 


summation  extend  over  finitely  many  terms. 

If  the  sequence  {m"}  converges  to  a  finite  limit  m  and  if  converges 

in  distribution,  then  Theorem  3.2  asserts  that  the  sequence  of  <!>' -valued  pro¬ 
cesses  converges  weakly  to  an  Omstein-Uhlenback  process  C.  with  diffusioi 

2 

rate  o  (given  by  (4.6)),  relaxation  rate  L  and  an  additional  drift  coeffi¬ 
cient  m. 

In  particular,  assuming  p^,  q^  to  be  finite  and  taking  m’’  =  m"  we  obtain, 
in  a  slightly  more  general  form,  the  main  result  of  [l4]. 

Example  2.  This  example  includes  the  cases  treated  by  Wan  and  Tuckwell  and 
by  Walsh  [21,  20].  Let  X  be  the  interval  [0,b].  If  X  represents  the  membrane 
of  a  neuron,  it  is  natural  to  consider  the  contraction  semigroup  {T^}  whose 
generator  -L  satisfies 


{/jl  s 


for  smooth  ^  with  compact  support  in  (0,b);  6i0  represents  the  rate  of  ion 
diffusion  within  the  neuron,  020  the  rate  at  which  ions  leak  across  the  morabrane 
In  the  cable  equation  approximation  to  the  electrical  properties  of  neuronal 


membranes 


1  2 


(see  [  7  ],  for  example),  5  =  ‘ 


6  =  ^  sec  ^  where  and  Rj  are  the  external  and  internal  longitudinal 

resistances  (in  ohm/cm) ,  R^^  the  membrane  resistance  per  unit  length  (in  ohm-cm) 

and  C,,  the  membrane  capacitance  per  unit  length  (in  F/cm) .  If  we  select  Le- 

besgue  measure  for  F  and  impose  Neumann  boundary  conditions  at  0  and  b  (i.c., 

seal  and  insulate  the  ends)  then  L  and  {T  }are  each  uniquely  determined  and 

2 

self-adjoint  on  H  =  L  (X,r) .  L  has  eigenfunctions  {(|)j}  with  associated  eigen¬ 
values  {X,}  given  by 


4»j(x)  =  Cj  cos(Trjx/b) 


j=0,l,2.. . 


If  we  set  C 
set  in  H. 


=  B  +  6(Trj/b)^. 

Assumption  A1  is  satisfied  for  any  >  ki  for  example,  yields 


Q  =  b  and  *  2\  for  j  2  1  then  {(|>j}  is  a  complete  orthonormal 


0.  =  =  ’sb  ((B+Dfi)"**  coth(b(B+l)VS  +  >5(B+1)"^ 

j  ^ 

By  Lemma  4.1  (with  r=0) ,  A2  is  satisfied  by  any  measure  u  onRx[0,bl  satisfying 
(4.2);  for  sufficiently  smooth  m,  we  may  take  in  (2.3).  With  the  choice 

of  X,  H  and  as  above  consider  the  Wan-Tuckwell  model  which  assumes  white 
noise  current  injection  at  a  single  point  x^  e  [0,b],  [  2ll»  The  impulse 
arrival-rate  measures  p”  are  taken  to  be  of  the  form 


(4.9)  /(AxB)  =  [  I  f^’" 

k*—  1  1 


(4.9)  is  a  special  case  of  (2.8).  Here  and  throughout  this  example  we  will 
assume  p,  q  <  “»  and  independent  of  n.  We  have 


where  and  ip  e  9  and 


(4.10)  =  ?  fg’"(a^’")^  > 

k=l  ®  ® 

Under  the  assumption 
(4. 11a)  lim 

IfW 

and 


Ji.n,  A,n.2 

^7  (a/)  • 


1  '  1 


(0  s  a  <  “) , 


(4.11b)  lim  max{a*^’",  af'”}  =  0, 
irw  k,Jl  ®  ^ 


a 


4.6 


the  assumptions  (A4)  (iii)-(v)  are  easily  seen  to  be  satisfied  with 


(4.12) 


The  white  noise  processes  then  have  the  representation 

(4.13)  x"(<|.)  =  t[m"((|,)  -  S"(<),)]  +  f  a^*"  /  (l>(x)N^'"(d.s.dx) 

k=l  ®  Xx{0,t]  ® 

-  ?  4’"  /  Nj’"(ds,dx), 

£=1  ^  X>‘(0,t]  ^ 

where  N^'”,  are  independent  Poisson  measures  with  variance  measures 

respectively  given  by  f^*"*v(dx)  and  f!’’”  v(dx)  with  v(B)  =  l-(x^).  In 

0  1  D  0 

(4.13)  in*'  is  given  by 


(4.14) 

with 


S"(<l>)  = 


**  k=l  ®  ®  £=1  ^  ^ 


Consistent  with  the  Wan-Tuckwell  model,  choose  m**(*)  =  ni*'(*)  and  assume  that 


(4.15) 


lim  Y  =  Y»  Y<“  exists. 
irx»  " 


-n 


Furthermore,  assume  5q  to  satisfy  A4  (i) .  Under  the  assumptions  (4. 9) -(4. 11) 
and  (4.15)  it  follows  from  Theorem  3.2  that  the  processes  satisfying 


dC"  =  -L'Cj  dt  +  dx",  (t>0) 


with  Cq  as  above,  converge  weakly  to  the  -valued  Gaussian  process  given 
by 


't  ■  ^ 


C4.16) 


The  real-valued  process  satisfy  the  stochastic  differential  equation 


(4.17)  dcj  =  +  Y<l>j(xQ)ldt  +  a|(j)j(Xp)  |dw^ 

with  initial  value  Each  is  a  standard,  real-valued  Wiener  process. 

Equation  (4.17)  is  equation  (25  )  obtained  by  Wan  and  Tuckwell  in  [21  ].  We 
note,  however,  that  in  [  21  ]  p  and  q  are  both  taken  to  be  equal  to  1.  As  a 
consequence  of  Theorem  2.3  it  also  follows  that  the  series  in  (4.16)  converges 
uniformly  in  0  s  t  ^  T  in  the  H  topology  where  q  >  max(^,-r  ). 

In  their  paper  Wan  and  Tuckwell  briefly  comment  on  an  extension  of  their 
model  in  which  the  impulses  can  arrive  at  different  points  of  X,  i.e.,  instead 

of  all  the  impulses  arriving  at  a  single  point  x^,  there  are  distinct  points 

k  H 

x^  (k=i,...,p)  at  which  excitatory  impulses  arrive  and  distinct  points  x^ 

(il=l,...,q)  which  receive  only  inhibitory  impulses.  The  authors  state  that 

the  locations  of  the  excitatory  and  inhibitory  synapses  "do  not  differ  by 

very  much". 

In  place  of  (4.9)  we  now  must  take 


(4.9)'  u"(AxB)  «  I  f^’"  *  ?  ^i'" 

k=l  Jl=l 


rhen 


k=l  *  ®  ®  ®  Jl=l  ^  ^ 


and  the  corresponding  m*'(^)  is  of  the  form 

,,  -n,..  f  k,n  k,  9  ^,n  l,n  ..  1. 

(4.18)  m  (♦)  =  i  fg  a^’  <J)(Xg)  -  i  a^’  (|)(x^). 


k=l 


4=1 


As  before,  we  will  take  m"  =  in”  in  the  choice  of  the  noise  process  x”.  The 
assumption 


Q  (♦»’!')  Q(<l>.4')  for  all  in  $ 


n-*<» 

,  ,q.  We  then  have 


nr+oo 


leads  to  the  existence  of  the  finite  limits  (assumed  positive) 
(4.19)  =  lim  f^’"(a^*")2,  =  lim  f?'»"(a^'")^  for  k=l,. 

and  i=l, . . 

=  f  4>(x^)'<'(x^)  +  ?  af  <(>(xf)ij/(xj'). 

1 _ lO  e  e  1  1 


k=l 


£=1 


This  assumption  and  the  additional  condition  (4.11b)  verifies  (A4)  (iii)-(v). 
(A4)  (i)  takes  the  form 


(4.20)  lim  ni”((|>)  =  m((j)),  finite  for  every  (J>  e  4>. 

n-»oo 

However,  conditions  (4.19),  (4.11b)  and  (4.20)  are  incompatible  unless 

(4.21)  P  ®  q  and  x*^  *  x!^  (k=l,...,p). 

G  i 

00 

(4.21)  can  be  easily  shown  by  using  the  fact  that  C  functions  with  compact 

support  belong  to  $.  Choose  such  a  function  with  'l>(x^)  =  1,  *l>(x™)  =  ((•(x^)  =  () 

1^ 

for  all  Z  «  l,...,q  and  mj^k,  where  x  is  a  point  not  belonging  to  the  set 

g  Ic  Ti  k  n 

{x.,  1=1,  —  ,q}.  Then  (4.20)  implies  f  ’  a  ’  ->  a  finite  limit  as  n  ->■  «> 

1  0  0 

k 

which  is  impossible  in  view  of  (4.19)  and  (4.11b).  Thus  all  the  points  x 

a  % 

belong  to  the  set  {x^^,  1=1,..., q}.  Similarly,  all  the  x.  belong  to  the  set 

{x^,  k=l,...,p}  and  thus  the  two  sets  are  the  same.  Renumbering  the  points 

i  k  k 

{x.}  if  necessary,  we  have  (4.21).  Writing  x  =  x.  =  x.  (k=l,...,p)  we  have 

X  ©IK 

Q(<fr,4>)  =  +  cx^) 

and 

m(<J>)  =  f  <l>(x. ) 


for  all  <p  £  9  where  y  =  lini[r  a  -  f .  a.  ],  Note  that  the  argument 

n-w  ^ 

just  given  holds  for  any  X  and  which  contains  functions  that  distinguish 
points  of  X.  The  stochastic  differential  equation  for  the  processes  corre¬ 
sponding  to  (4. 17)  now  takes  the  form 

(4.22)  d^J  =  +  i(<|.j)]dt  ♦  /Q((|)j.(j).)  dwj 

with  initial  value  and  w^  a  standard  Wiener  process.  The  4>' -valued  process  C 
is  given  by  (4.16)  with  satisfying  (4.22).  The  processes  5*^  satisfy  an  equat 
tion  very  similar  to  that  obtained  in  the  simpler  Wan-Tuckwell  model  with  white 
noise  current  injection  at  a  single  point.  Equation  (4.22)  differs  from  (4.17) 
only  in  the  constants  m(<|)^)  and  i/Q(<j)^  ,(j)^)  . 

Thus  if  C*  is  a  solution  to  the  initial  value  problem 


dC^  =  -L'C^dt  ♦  dX^ 


^0  *  ^0 


in  which  is  a  generalized  Poisson  process  with  only  very  small  jumps,  then 
C,  is  close  in  distribution  to  the  solution  C,  to  the  initial  value  problem 


dc^  =  -L'^^dt  +  dW^ 


^0  "  ^0 


If  W^  is  a  Gaussian  ^'-valued  process  with  the  same  mean  and  covariance  as  X^ 
and  if  Cq  is  a  Gaussian  -valued  random  variable  whose  distribution  is  close 
to  that  of  This  supports  the  use  of  Gaussian  methods  to  study  the  approx¬ 
imate  distribution  of  the  process  including  first-passage  times  and  re¬ 
lated  functionals  (e.g.  in  [21]  and  [ig])- 
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In  [20  ]  Walsh  considers  a  model  in  which  X,  F  and  {T^}  are  the  same  as 
in  the  Wan-Tuckwell  example  discussed  above  but  the  choice  of  y  is  more  gen¬ 
eral,  i.e.,  y  is  determined  by  a  finite  measure  v  on  X  and  a  Markov  kernel 
K  via  the  relation 

yCAxB)  =  /  K(x,A)  v(dx),  A  c  K,  B  c  X. 

B 

The  kernel  K(x,da)  is  the  regular  conditional  probability  distribution  of  the 

size  a  €  R  of  impulses  arriving  at  the  site  x  e  X,  while  the  measure  v(B) 

is  the  overall  arrival  rate  yQRxB)  of  impulses  of  all  sizes  at  points  xc  B  c  X 

2  2 

The  requirement  made  in  [20  ]  that  /[/a  K(x,da)  +  (fa  K(x,da))  ]  v(dx)  < 

entails  that  4.2  is  satisfied  so  that  A2  holds  and  Theorems  2.1  and  2.2  apply. 

Incidentally,  Walsh  represents  T^  through  its  integral  kernel  G(x,y;t)  = 

-tX. 

e  ^  <l>j  (x)  <l>j  (y) .  The  series  converges  uniformly  in  x  and  y  by  Al.  Walsh 
defines  a  two-parameter  stochastic  process  V(t,x)  which  is  an  integral  kernel 
for  our  : 

=  /  <l>(x)  V(t,x)  r(dx) 

and  our  notation  is  related  to  his  as  follows: 

VgCxjrcdx), 

X 

=  /  <l>(y)  V(t.y)r(dy) 

^  X 

=  /  <l>(y)  /  v»(x)  G(x,y;t)r(dx)r(dy) 

X  X 

t 

+  /  <P(y)  [f  f  G(x,y;t-s)  F(dsdx)  ]r(dy) 

X  X  0 

=  /[/  G(x,y;t)(f)(y)dy]  v  (x)r(dx) 

X  ^ 

t 

*  /  /[/  G(x,y;t-s)({i(y)r(dy)]  F(dsdx) 

0  X 


=  /[V(x)]  VQ(x)r(dx) 
X 

t 


W  /  [T 

0  IRxX 


t-s 


4>(x)]  aN(dadxds) 


(O.t]  ^  ^ 


It  is  shown  in  Theorem  3.3  of  [20  ]  that,  under  certain  conditions 
(3.1  of  [20]),  almost  all  of  the  finite- dimensional  distributions  of  a 
sequence  of  such  processes  v”  converge  to  those  of  a  Gaussian  process  V  which 
is  ai\  integral  kernel  for  our  Walsh's  condition  (3.1a)  implies  our  A4  (iii) 
while  his  (3.1b)  is  equivalent  to  our  A4  (iv) ' .  It  should  be  noted,  though, 
that  our  approach  to  the  problem  is  different  from  Walsh's.  Our  results  are 
not  concerned  with  the  convergence  in  distribution  of  v”(t,x)  for  fixed  t 
and  x.  On  the  other  hand.  Theorem  3.2  of  our  paper  stengthens  Theorem  3.3 
of  [20]  by  proving  the  convergence  in  distribution  of  the  $' -valued  processes 
C"  in  the  sense  that  the  probability  measures  in  D(1R^:'I>')  induced  by  the  ^"'s 
converge  weakly  to  that  induced  by  The  advantage  of  establishing  weak 

convergence  of  the  -valued  processes  is  that  it  yields  convergence  in  distri¬ 
bution  of  continuous  functionals  of  the  paths,  such  as  the  first  passage  times 
t(4))  =  inf{t>0:  C^(<j>)>0},  =  “>  if  {...}  is  the  null  event,  ((|)  e  $  or  H^) . 

As  already  noted  earlier,  our  approach  to  modelling  the  biological  phen¬ 
omena  of  interest  as  -valued  processes  is  very  general  as  is  clear  from  the 
assumptions  made  in  Theorem  3.1  on  X,  F  and  the  semigroup  {T^},  and  as  will 
be  further  illustrated  in  the  next  example  and  Section  5.  In  these  examples, 

X  will  be  a  sphere  in  R  or  a  smootlv  compact  manifold  with  or  without  boundary 


.Such  problems  lie  outside  the  scope  of  multi-parameter  martingale  theory.  Even 


2 


when  X  cR  (d^l),  the  semigroup  {T^}  cannot,  in  general,  be  represented  l)y 
a  kernel  given  by  a  uniformly  convergent  series,  the  case  treated  by  Walsh. 
When  the  series  makes  sense  as  a  distribution,  one  is  naturally  led  to  a 
generalized  stochastic  process  of  the  kind  we  have  considered  in  this  paper. 


Example  3 .  X  consists  of  a  sphere. 

2  3 

Let  X  =  S  ,  the  unit  sphere  in  R  ,  with  Lebesgue  surface  measure  P.  If 
Ag  denotes  the  spherical  Laplace-Beltrami  operator 


2 

(where  0,ri  are  the  Euler  angles  on  S  ),  then  L  =  -6  +  6A  is  self-adjoint 

B 

2 

on  H  =  L  (X,r)  for  any  real  numbers  B,  5.  If  B  and  6  are  positive  then  -L 
generates  a  contraction  semigroup  {T^}  as  before.  This  time  the  eigenfunctions 
are  the  spherical  harmonics  (£=0,1,...;  m=-£,...,£)  with  eigenvalues 
A  =  B+6£(£+l)  for  L,  e  for  T^;  we  set 

4>j  =  j=m+Jl(Jl+l)  (i.e.  £=[Ij^*]]  ,  m=j-£^-£) 

A.  =  6+6  aj**]]  ([[j’^D+D 


Assumption  A1  is  satisfied  for  any  r^  >  *5,  since 

-2r 


-2r, 


^  =  J;^(2£+1)(1+B+6  (£+£))  ^ 


-2r  -2r 

<  4B  +6  V(r^-»S) 


<  «. 


The  spherical  harmonics  obey  the  bound 
sup  lY,  (x)l^  =  (2£+l)/4Tr 


so  {<!> 


.}  satisfies 
J 

sup  sup  l<l>.(x)|  (1+X.)  ^  =  (4iT^min(6,43))  ^ 

j  X  ^  ^ 


By  Lemma  4.1,  any  y  on  R><X  satisfying  (4.2)  also  satisfies  A2.  For  further 
details  the  reader  may  consult  [5  ]. 


Applications  Where  X  is  a  Compact  Riemannian  Manifold 


Often  the  set  X  is  intended  to  be  the  cell  wall  or  surface  membrane  of  a 
neuron;  as  such  it  ought  to  admit  a  mathematical  representation  as  a  compact 
Riemannian  manifold  of  two  dimensions.  In  voltage-clamp  experiments  the  cell 
wall  is  cut  and  the  membrane  potential  voltage  at  the  cut  is  held  at  some  nre- 
scribed  value  by  the  experimental  apparatus;  in  this  case  X  may  be  regarded  as 
a  manifold  with  boundary. 

The  following  theorem  gives  conditions  on  a  manifold-with-boundary  X,  a 

2 

measure  T  on  X,  and  a  semigroup  {T^}  on  L  (X.T)  which  guarantee  that  assump¬ 
tions  Al)  and  A2)  will  be  satisfied  for  any  measure  p  satisfying  (4.2),  i.e. 
for  any  noise  process  admitting  mean  and  covariance  measures.  See  Hbrmander 
fril  ),  ch.  X]  for  details  about  elliptic  boundary  systems  on  manifolds. 

Theorem  5.1 

Let  X  be  a  smooth  d-dimensional  compact  Riemannian  manifold  with  smooth 

(possibly  empty)  boundary  3X  and  Riemannian  volume  element  dP.  Let  L  be  a 

2 

positive  self-adjoint  operator  on  a  domain  PcH=L  (X,r)  satisfying 

i)  c”(X)cl7 

00 

ii)  The  restriction  L_  of  L  to  C  (X)  is  a  uniformly  strongly  elliptic 

V  ^ 

differential  operator  of  order  2m >  0  with  smooth  coefficients. 

iii)  Hilbert  space  of  those  elements  in  H  with  2m  weak 

derivatives  in  H. 

Then  L  admits  a  complete  orthonorraal  set  (4)^}  of  eigenfunctions  in  H  with 
eigenvalues  {X^}  satisfying 
i )  L(J) .  =  A  .  4) . 

J  r  .1 

ii)  4>j  e  C"(X) 


”  1 

iii)  Ij(l+^j)  ^  “  for  all  rj^  >  d/4m 


iv)  sup  .sup  l4>^(x)  1  (1  +  X.)' 
j  X  ^  ^ 


V  >  •  •  ■ 

'•  A  • 


<  00 


for  all  r  >  d/4m. 


Corollary  5.1 

The  contraction  semigroup  generated  by  -L  satisfies  Al.  Any  measure  v* 
on  X X  satisfying  (4.2)  must  also  satisfy  A2. 

Proof  of  Theorem  5.1 

In  case  X  is  an  open  set  in  ,  Theorem  14.6  of  Agroon  ri]  gives  i)  and 
iii);  ii)  follows  from  elliptic  regularity  (e.g.  ril  Theorem  9.3)  and  iv)  from 
the  .Sobolev  imbedding  theorem  (e.g.  fl]  Theorem  3.9).  Completeness  of  {<}>^} 
follows  from  r2],  Theo7*em  3.4. 

The  case  of  a  general  Riemannian  manifold  can  be  treated  by  extending 
Agmon's  proof  to  manifolds  (using  the  Schauder  estimates  of  <‘31)  or  by  employing 
a  partition-of-unity  argument  as  in  Hormander  (Til)  Chapter  X). 

Remark  5.1 

Frequently  in  applications  Lq  is  given  as  part  of  an  elliptic  boundary  sys¬ 
tem  L=  {1,^;  (see  ’’ll],  p.  254)  in  which  are  smooth  dif- 

00 

fercntial  operators  of  order  <  2m  on  C  (3X) .  Theorem  5.1  will  apply  if  it  can 

go 

be  shown  that  1,^  is  essentially  self-adjoint  on  the  domain  {(}>e  C  (X): 
has  a  continuous  extension  ((>  to  a  =  Xu3X,  the  restriction  of  (j>  to  3X  is  C  ,  and 
^j(f)~U  on  3X  for  l<j<m}. 

The  second-order  elliptic  operators  with  Dirichlet  or  Neumann  boundary 

conditions  can  be  treated  in  this  way. 

examples  1,  2  and  3  of  Section  4  are  special  cases  of  this,  with  d  =  0, 

(1-1,  and  d  =  2  respectively.  A  more  realistic  example  would  appear  to  be  that 

in  which  X  is  the  boundary  of  a  3-dimensional  solid,  so  9X  =  0,  and  L  = -A„  + 

B 

3(x)  is  the  sum  of  the  Laplace-Beltrami  operator  and  a  smooth  positive  function 
representing  the  trans-membrane  conductance  on  a  (not  necessarily  spherical)  2- 

li 


manifold  X, 


Example  4. 


Before  simplifying  to  X=ro,b]  Wan  and  Tuckwell  discuss  a  one-dimensional 
finitely  branching  tree  as  a  model  for  X: 


Although  this  is  not  quite  a  manifold  (it  is  not  locally  homeomorphic  to  K  at 
the  nodes)  it  nonetheless  can  be  handled  by  our  methods  as  follows. 

Represent  X  as  the  disjoint  union  of  N  compact  intervals  ro.bji . 

and  each  function  4i  on  X  as  a  family  of  functions  on  the  intervals 

satisfying  certain  boundary  conditions.  If  represents  the  (possibly  empty) 
set  of  indices  for  the  segments  emerging  from  the  node  at  the  end  of  the  ith 
segment ,  the  appropriate  boundary  conditions  are 


1)  <t>^(b^)  =<)>j(0)  for  each  je  IsisN 

2)  (t).'(b.)=  I  (|):(0)  for  each  l<isN  , 

^  ^  -5.-  T  J 


3)  <|.|(0)  =0  . 

The  first  of  these  guarantees  that  <p  will  be  continuous  when  regarded  as  a 
function  on  the  tree,  while  the  second  is  an  expression  of  Kirchhoff’s  I, aw 
forbidding  any  leak  of  current  at  the  nodes.  The  third  equation  imposes  Neumann 
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(or  "insulating")  boundary  conditions  at  the  base  of  the  tree;  it  or  any  of 
the  N  conditions  in  2)  could  be  replaced  by  a  Dirichlet  (or  "grounded")  condi¬ 
tion.  It  is  easy  to  integrate  by  parts  and  show  that  A(|)^  =  (|)V  is  self-adjoint 
on  the  domain  P  of  all  functions  on  X  satisfying  the  boundary  conditions 
above  and  possessing  an  absolutely  continuous  derivative  and  a  square-integra- 
ble  second  derivative  with  respect  to  Lebesgue  measure  F. 

Theorem  5.1  still  doesn't  quite  apply  since  the  boundary  operators  are 

non-local  (on  the  disconnected  set  X),  but  its  conclusion  still  holds  true. 

00 

In  fact,  -A  + 6  has  a  uniformly  bounded  complete  orthonormal  set  of  C  eigen¬ 
vectors  and  hence  one  can  take  r  =  0  instead  of  just  r  s 1  in  conclusion  iv) . 


Nonlinear  Problems  and  Directions  for  Future  Work 


The  stochastic  behavior  of  the  voltage  potential  investigated  in  the  previous 
sections  has  the  form  of  a  differential  equation 

(6.1)  dC^  =  -L'?^dt  +  dX^ 

where  L'  :  is  the  adjoint  of  L.  In  our  model  L'  is  a  linear  operator  and 

hence  {Tp  appearing  in  (3.17)  is  a  linear  evolution  semigroup.  This  is 

a  decidedly  unrealistic  feature,  for  physically  it  amounts  to  the  assumption  that 
the  electrical  properties  of  the  membrane  are  unaffected  by  changes  in  the  poten¬ 
tial  voltage  across  the  membrane.  Our  model  fails  to  take  into  account  experi¬ 
mentally  observed  features  of  certain  neurons,  for  instance  the  following: 

(a)  As  the  cell  membrane  is  progressively  depolarized,  the  postsynaptic  poten¬ 
tials  (inhibitory  as  well  as  excitatory)  eventually  become  reversed  in  sign. 

The  reversal  potential  would  introduce  one  type  of  nonlinearity  in  the  model  (sec 

[21]) 

(b)  At  least  in  some  cells  (such  as  sympathetic  ganglion  cells,  see  [4], 
p.  135)  the  sizes  of  the  postsynaptic  potential  impulses  are  dependent  on  the 
state  of  depolarization  of  the  membrane  potential. 

It  would  thus  seem  that  a  model  which  better  describes  the  physiological 
process  would  let  satisfy  an  equation  such  as 

(6.2)  dC^  =  a^dt  +  b^dX^ 

in  which  a'  is  a  ^’-valued  process  adapted  to  {X^}  (most  likely  a  nonrandom, 
nonlinear  ^'-valued  function  of  5^)  and  b^  is  an  adapted  process  taking  values 
in  a  space  of  functions  from  4>'  to 

We  can  replace  X^  by  W,  and  K,  by  C,  in  (6.1)  and  (6.2)  and  seek  to  approxi¬ 
mate  the  behavior  of  the  process  by  that  of  a  nuclear  space-valued  diffusion 
satisfying 

(6.3)  dC.  =  aMt  +  bMW^  , 


6.2 
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but  the  techniques  used  in  the  present  paper  do  not  suffice  to  prove  the  existence 
and  uniqueness  of  solutions  to  (6.2)  and  (6.3)  or  to  prove  the  distributional  con¬ 
vergence  of  a  sequence  of  solutions  to  (6.2)  to  the  corresponding  solution  of 
(6.3)  when  a'  and  b'  are  nonlinear. 

Even  the  use  of  a  semigroup  (linear  or  not)  {T^}  to  model  the  evolution  of 
the  membrane  potential  in  the  absence  of  incident  impulses  entails  an  implicit 
assumption  about  the  physical  system  which  is  unrealistic,  namely  that  the  state 
of  the  system  is  unambiguously  specified  by  giving  only  the  membrane  potential 
itself  at  every  point  of  X.  It  is  known  [13]  that  the  local  membrane  behavior 
depends  critically  on  the  concentration  gradients  of  sodium,  potassium,  chloride, 
and  calcium  ions,  and  that  active  transport  of  these  ions  as  well  as  diffusion 
play  important  roles.  Recent  experiments  revealing  the  stochastic  behavior  of 
ion-specific  gates  through  cell  membranes  offer  new  opportunities  for  more 
elaborate  stochastic  modelling  of  this  important  system. 

Finally,  the  stochastic  differential  equation  model  with  which  this  paper 
deals  and  its  extensions  briefly  indicated  above  have  applications  to  other  areas 
of  biology,  e.g.  to  problems  of  emigration  of  biological  populations.  We  hope  to 
investigate  such  problems  in  future  papers. 
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